yAD-A036  523 


unclassified 


CLEMSON  UN IV  S C COLL  OF  ENGINEERING  F/G  13/2 

A SYSTEMS  ANALYSIS  OF  WATER  QUALITY  SURVEY  DESIGN.  APPENDIX  II.— ETC(U) 
AUG  75  L C WILCOX.  B E GILLILAND  DADA17-72-C-2152 


I*  2 

A0A038523 


A SYSTEMS  ANALYSIS  OF  JfATER  QUALITY  .SURVEY RESIGN 


APPENDIX  JXT- 
DOCUMENTATION 


§URVEY_PLANNTN6  ^PROGRAM  user  s manual 


SUPP0RTEDBY 

ARMY  MEDICAL  RESEARCH  AND  I 
WASHINGTON-,  D.C,  , 


.OPMENT  COMMAND 


PUBLICATION  BY 

ENGINEERING  RESEARCH 
CLEMSON  UNIVERSITY 
CLEMSON,  S,C. 


A SYSTEMS  ANALYSIS  OF  WATER  QUALITY  SURVEY  DESIGN 


Documentation  of  the  Survey  Planning 
Computer  Program.  A User's  Manual 


Authors:  Dr.  L.  C.  Wilcox,  Dr.  B.  E.  Gilliland, 
Dr.  R.  W.  Gilchrist,  Dr.  T.  M.  Keinath 


August  1975 


Supported  by 

U.  S.  Army  Medical  Research  and  Development  Command 
Washington,  D.  C.  203 1 ^ 

Contract  No.  DA  DA  1 7-72-C-2 1 52 

This  document  has  been  approved  for  public  release; 
distribution  unlimited. 


The  findings  of  this  report  are  not  to  be  construed 
as  an  official  Department  of  the  Army  position  unless 
so  designated  by  other  authorized  documents. 


D D C 


f*P  »r  /'«■/• 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1 Ni.r-OH  T NUMFiKH 


Final  Report 
*.  title  (and  Subtitle) 


\2.  GC  r AC  F S N K&.I  3 WrCPlLsT’S  C*T*i  OG  t.  :i  >J  R 


5.  TYPE  OF  REPORT  4 P.  R|OD  COVERED 


A Systems  Analysis  of  Water  Quality  Survey  Design 


7.  author^;  1 

Lyle  C.  Wilcox  / Bobby  E.  Gilliland 

Thomas  L.  Drak£  / 

R a 1 p h W.  Gil Christ  ' 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Clemsorl  University 
College  of  Engineering 
Clemson,  SC  29631 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Final 

6.  PERFORMING  ORG.  REPORT  NUMBER 

Final 

8.  CONTRACT  OR  grant  NUMBERf*) 

DADA  17-72-C-2152  . 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4 WORK  UNIT  NUMBERS 


12.  REPORT  DATE 


August  1975 
"TT  number  of  p ages 


U MONITORING  AGENCY  NAME  4 ADDRESSf/T  dltterenl  from  Controlling  Ottlce)  15.  SECURITY  CLASS,  (ot  thl,  report) 


IS*.  DECLASSIFICATI  ON/ DOWNGRADING 
SCHEDULE 


16-  DISTRIBUTION  STATEMENT  (a!  I hit  Report) 


Distribution  of  this  document  is  unlimited. 


'7.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  ftlock  20,  It  different  from  Report) 


18  SUPPL  i MENT  ARY  NOTES 


19.  KEY  WORDS  (Conti. :ue  on  reverse  side  If  necessary  and  Identify  by  block  number) 


Systems  Analysis 
Water  Quality  Modeling 

Water  Quality  Survey  Simulation 

Water  Quality  Survey  Design  Ammunition  Manufacturing  Processes 

• — A rmy -Ammuft4t+en-P-lant^ W^-r  Quail tyjteasuremen t s 

20.  AL'iTHACT  (Continue  o n.  reveree  s,da  ft  neceaenry  and  Identify  by  block  number) 

This  is  the  final  report  of  a three  year  project  titled,  "A  Systems  Analysis 
of  Water  Quality  Survey  Design." 

In  this  project  a study  was  made  of  water  quality  surveys  conducted  hy  the 
United  States  Army  F.nvi  ronmental  Hygiene  Agency  (AEHA).  Mainly  data  and  reports 
from  studies  of  Army  Ammunition  Plants  (AAP)  were  used. 

The  focus  of  this  project  was  the  deve1opment  of  computer  aided  procedures 
which  would  assure  efficient  use  of  manpower  and  equipment  and  assure  that  the 
measurements  taken  eive  a reasonalbe  representation  of  the  system.  Planning 


DD  I 1473  EDITION  OF  > NOV  «S  IS  OBSOLETE 

J Unlimited 

FFCURITY  CLASSIFICATION  of  Tm’S  PAC.F  (KTttn  D«f«  Fntereef) 


LCUKlTY  CLASSIFICATION  Of  THIS  PAGEf'*>«n  pati  I rtrrad) 


survey,  conducting  the  survey  and  reporting  on  the  survey  were  included  in  the 
study . 

The  site  modeling  program  models  the  manufacturing  processes  which  contri- 
bute pollutants  to  the  system,  models  the  sewer  system,  and  models  the  treat- 
ment system  including  acid  or  caustic  neutralization,  settling  ponds,  and 
domestic  treatment.  The  inputs  to  the  model  are  the  production  levels  of  the 
manufacturing  processes  and  the  outputs  are  the  predicted  pollutant  measuremen 
values  at  each  possible  measure  point  in  the  system. 

The  resource  matching  program  accepts  data  defining  proposed  measurements 
and  matches  these  against  the  available  time,  manpower,  and  equipment.  The 
output  lists  the  pollutant  to  be  measured  at  each  measure  point,  the  total 
commitment  of  time  for  each  analyst  and  for  each  piece  of  equipment.  Note  is 
made  of  any  overcommitment  of  manpower  or  equipment. 

The  model  refinement  or  updating  program  accepts  measurements  taken  durinj 
a preliminary  survey  or  during  a regular  survey  and  computes  suggested  new 
parameters  for  the  process  models. 

\ The  indicator  model  program  evaluates  the  performance  of  sanitary  treat- 
men  tVfacili  ties . 

The  program  uses  design  data,  data  from  the  operating  log  and/or  data 
generated  during  the  survey  and  computes  key  operational  characteristics.  Com- 
paring these  with  desirable  values  as  cited  in  design  books  and  manuals  will 
give  the  survey  planner  insight  into  the  operation  of  the  system  and  suggest 
the  need  for  more  survey  measurements  or  the  need  for  changes  in  operation. 

A system  was  developed  for  automatic  instrumentation  of  pH,  conductivity, 
and  other  parameters  which  use  strip  charp  recordings.  Interface  hardware  was 
selected  and  purchased  and  interface  software  was  developed  for  direct  connec- 
tion to  a digital  computer. 

A data  handling  system  was  developed  for  use  during  and  after  the  survey. 
A PDP8-OS/8  and  peripheral  equipment  was  purchased.  Software  was  developed  to 
perform  data  handling  functions  and  to  direct  the  user  in  application  of  the 
software.  The  program  accepts  raw  data  from  the  analytical  chemist  and  per- 
forms data  conversions,  transcriptions,  and  data  logging  functions.  Output 
is  available  in  several  forms  as  may  be  needed  for  various  reports  during  and 
at  the  end  of  the  survey.  i 

Recommendations  are?  the  survey  planner  should  obtain  sufficient  data  in 
a preliminary  survey  to  model  and  analyse  the  site;  measurements  should  he  auti 
mated  to  the  maximum  extent  possible;  data  handling  should  be  delegated  to  tin- 
computer  when  the  operations  are  well  defined  and  repetitive.  The  programs, 
software  and  hardware  included  here  will  assist  the  survey  planner  in  followin 
these  recommendations  and  design  a more  effective  survey. 
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INTRODUCTION 

This  manual  was  written  to  assist  in  the  use  of  the  Survey  Planning 
Program  developed  as  part  of  the  research  entitled  "A  Systems  Analysis  of 
Water  Quality  Survey  Design"  conducted  by  the  Engineering  College  at  Clem- 
son  University  and  supported  by  U.S.  Army  Medical  Research  and  Development 
Command . 

The  program  can  perform  three  basic  functions  by  assembling  the  proper 
subroutines.  These  functions  are: 

1.  Modeling  and  simulation  of  the  plant  site  operation.  The  model 
includes  the  effluent  of  the  manufacturing  processes,  the  sewer 
system  pollutants,  and  industrial  and  domestic  sewe^;  treatment. 

2.  Resource  matching  or  resource  commitment  for  a planned  survey. 
Using  the  modeling  and  simulation  results,  the  survey  planner 
determines  the  frequency  and  location  of  the  measure  points  for 
each  pollutant  to  be  measured.  The  RM  model  uses  this  informa- 
tion along  with  information  on  available  resources  (time,  per- 
sonnel and  equipment)  determines  the  required  commitment  of 
personnel  and  equipment.  Note  is  made  of  non-feasible  survey 
plans  whioii  would  require  more  personnel  or  equipment  than  is 
available.  The  survey  planner  can  then  add  more  oersonnel  or 
equipment  for  the  survey  or  try  another  survey  plan  which  will 
use  different  personnel  and  equipment.  By  iterating  through 
several  survey  plans  the  planner  can  find  a survey  plan  which 
will  yield  the  maximum  of  useful  data  for  the  available  time, 


personnel  and  equipment. 
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3.  Model  Verification  and/or  Updating. 

The  site  model  process  models  can  be  updated  or  improved  using 
preliminary  survey  data  or  survey  results.  While  the  best  data 
available  were  used  in  constructing  the  process  models,  there 
was  considerable  variation  in  data  for  a given  process  from  one 
plant  to  another  or  from  one  season  to  another.  Better  process 


model  values  may  be  obtained  using  the  CORRCT  subroutine  with 
judiciously  selected  grab  samp  1 es.  taken  on  a preliminary  survey. 
Survey  results  can  also  be  used  to  verify  or  refine  the  process 
models.  With  confidence  in  the  model,  it  can  now  be  used  to 


predict  levels  of  potential  pollutants  for  production  levels 
other  than  those  which  occurred  during  the  survey.  This  should 
give  added  significance  to  the  survey  report. 

The  site  modeling  and  simu'ation  uses  MAIN  and  the  following  sub- 
routines: PROCES,  START,  TOP,  NEUTRA,  SETTLE,  DOME. 

The  resource  matching  application  uses  all  of  the  subroutines  re- 
quired for  site  modeling  plus  the  following  subroutines:  LEVEL,  FLAG, 

RM,  CONCK,  TO,  EQCHER,  VSCHEC,  CNCHEC,  ADD,  SET,  PIEQCH,  SUBT,  INFORM, 
SAMPLE,  PRPAR , USECT,  TICHAN. 


The  model  verification  and/or  updating  uses  all  of  the  subrouintes 
required  for  site  modeling  plus  the  subroutine  CORRCT. 

The  functions  of  the  subroutines  are  discussed  briefly  below. 

MAIN:  sets  the  system  dimensions  and  calls  appropriate  subroutines. 

PROCES:  This  program  reads  in  process  mnemonics  for  processes  present 

at  the  particular  site  being  studied.  Appropriate  process  data  is 
called  out  of  process  library.  Process  data  can  be  altered  to  better 
represent  a particular  site. 


1 


3 


Process  Library:  Contains  process  models. 

START:  Reads  in;  the  number  of  processes (sources) , the  number  of 

branches,  the  number  of  parameters  and  the  number  of  outfalls  in 
this  site  study.  Also  reads  in  the  numbers  of  the  particular  pollu- 
tants in  the  current  study  (NPLIST).  Reads  in  branch  names.  Branches 
are  named  for  convenient  interpretation  of  output  data.  Calls  CHK1 . 

CHK1 : Reads  in  the  number  of  non-competing  parameters  (used  in  RM) ; 
reads  in  the  names  of  all  of  the  pollutants  in  the  master  list,  for 
each  the  number  of  measurement  methods  available  and  the  sampling 
frequency.  This  program  also  prints  out  the  input  information  available 
at  this  point.  In  particular  this  includes  the  process  data  and  the 
topology  matrix. 

TOP:  This  program  essentially  does  the  following:  reads  in  treatment 

data;  takes  process  data  and  site  topology  data  and  calculates  the  flow 
and  pollutant  parameters  for  each  branch  of  the  system.  The  effects 
of  treatments  are  included  in  the  calculation;  calls  the  treatment  sub- 
routine NEUTRA,  SETTLE,  and  DOME;  calls  CHK2  to  output  results. 

NEUTRA:  This  subroutine  is  called  by  TOP  if  acid  or  caustic  neutrali- 

zation is  required  in  any  branch. 

SETTLE:  This  subroutine  is  called  by  TOP  is  a settling  pond  or  clari- 

fier is  required  in  any  branch. 

DOME:  This  subroutine  is  called  by  TOP  if  domestic  or  sanitary  treat- 

ment i s requ i red . 

CHK2:  This  program  essentially  writes  out  the  results  of  the  TOP  cal- 

culations, i.e.,  the  system  simulation. 

FLAG:  This  subroutine  is  not  used  in  system  simulation  but  in  system 

planning  with  RM.  However,  it  is  called  by  TOP  to  prepare  data  for  RM. 
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LEVEL  arranges  TOP  output  for  use  in  RM. 


FLAG  allows  survey  planner  to  specify  for  RM,  measure  points  in 


addition  to  those  selected  by  RM  logic. 


RM  reads  in  and  prints  out  resource  information,  allocates  resources 


and  calls  subroutines  when  needed  to  relieve  violations,  and  prints 


out  resource  allocation  information. 


CONCK  determines  feasible  measure  methods  and  allocates  resources  to 


flagged  points. 


TO  zeros  out  temporary  allocations  for  equipment  times. 


EQCHEC  determines  if  any  equipment  time  violations  occur. 


VSCHEC  checks  for  van  space  constraint  violations. 


CNCHEC  checks  for  analyst  time  constraint  violation. 


ADD  adds  current  allocations  of  equipment  time,  analysts'  times,  van 


space  and  costs. 


SET  stores  the  array  of  methods  used  at  a branch  in  a level  for  a given 


pa  ramete  r . 


PIEQCH  searches  past  assignments  to  attempt  to  relieve  constraint  viola- 


tions. 


SUBT  subtracts  current  allocations. 


INFORM  prints  out  branch  names  and  parameters  to  be  measured  at  that 


branch . 


SAMPLE  allocates  analyst  type  1 time  required  for  taking  samples  and 


measuring  flow. 


PRPAR  outputs  parameter  measurement  information. 


USECT  total  usage  equipment, 


TICHAN  converts  minutes  to  hours  and  minutes. 


■ 


CORRCT  is  used  for  model  verification  or  updating. 

Subsequent  sections  of  this  manual  define  the  data  required  by 
these  program  subroutines,  list  the  required  coding  and  ordering  of  data, 
and  illustrate  the  format  and  computer  card  layout  required  by  the  sub- 
routine. Also,  the  various  components  of  the  output  information  are  listed 
and  explained  in  the  last  section  of  this  manual. 

The  programs  have  been  written  so  that  only  the  subroutines  required 
for  the  particular  application  need  to  be  read  in.  Thus,  the  program  could 
be  visualized  as  three  separate  decks  with  some  of  the  subroutines  appearing 
in  more  than  one  deck.  Also,  the  complete  set  of  subroutines  can  be  read 
as  one  big  single  program.  Depending  on  the  computer  facilities  available, 
memory  space  available  and  turn-around  time  may  be  a problem  if  the  three 
smaller  decks  are  not  used. 

In  either  case  flags  have  t:o  be  read  in  MAIN  to  define  the  particular 
application  that  is  desired. 

NFLAG  =0  call  LEVEL  and  RM 

= 1 terminate  after  TOP  or  CORRCT 

EPSLON  =0  do  not  call  CORRCT 

=1  call  CORRCT 
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BRIEF  DESCRIPTION  OF  FORMAT  STATEMENTS 

The  purpose  of  the  FORTRAN  FORMAT  statements  is  to  tell  the  com- 
puter what  kind  of  information  is  to  be  inputted  and  where  to  find  that 
information  on  a data  card.  In  this  program,  only  three  different  types 
of  data  are  inputted:  integer  data  ('T1  FORMAT),  floating  point  data 

("F"  FORMAT)  and  alphanumeric  data  ("A"  FORMAT). 

Integer  Data 

Integer  formats  are  of  the  form  wlm,  where  w is  the  number  of  ; 

integer  data  fieios,  and  m is  the  length  of  each  integer  fieid.  For 
example,  315  !M  read  3 items  of  data  from  a card,  with  each  item  in 
a field  5 card  columns  wide.  It  should  be  mentioned  that  within  a field 
the  data  must  be  right  justified.  This  means  that  i f a 14  format  is 
used  to  read  in  the  number  6,  6 must  appear  in  the  4th  column  of  the 
field.  If  it  does  not,  the  trailing  blanks  are  read  as  zeros.  So,  if 
the  6 is  placed  in  the  3rd  column,  it  will  be  read  by  the  computer  as 
60 , or  in  the  2nd  column  it  will  be  read  as  600.  Care  must  be  taken 
that  integers  are  right  justified.  Many  apparent  programming  errors  may 
be  caused  by  incorrect  punching  of  data  cards. 

Exampl es : 

READ  100,  K,  L,  M 


100  FORMAT  (314) 


This  will  read  in  three  integer  values  for  the  variables  K,  L,  and  M. 
If  the  desired  values  for  those  variables  are  5.  1 06  and  79.  the  data  card 
would  appear  as  follows: 


Col umn 


Since  the  computer  will  only  look  at  the  first  12  columns  in  this 
example,  columns  1 3“80  could  be  filled  with  anything. 


Floating  Point  Data 

Floating  Point  Formats  are  of  the  Form  wFm.n,  where  w is  the  number  of 
floating  point  fields,  m is  the  total  width  of  each  field,  and  n is  the 
number  of  places  to  the  right  of  the  decimal  point.  The  restriction  on  n 
is  n 4 m+1 , since  one  place  is  reserved  for  the  decimal  point.  Floating 
point  allows  the  user  to  omit  the  decimal  point,  i f the  data  i f right 
justified.  Otherwise,  data  may  be  placed  anywhere  in  the  Field. 

Example : 

READ  101 , X 
101  FORMAT  (F10.3) 

If  the  value  desired  for  X is  107.3^6,  the  data  card  may  appear 


as 


6 


or  as 


2) 


The  presence  of  the  decimal  point  overrides  the  Format  specification 
In  the  first  case,  if  1073^6  started  in  column  4 instead  of  column  5 it 
would  be  read  in  as  1073.^60. 


Alphanumeric  Data 


Alphanumeric  formats  are  of  the  Form  wAn,  where  w is  the  number  of 
alphanumeric  fields,  and  n is  the  width  of  each  field.  Any  character  may 
be  read  in  as  alphanumeric  data. 

Example: 

READ  102,  NAME,  N,  J 
102  FORMAT  (2AA.A3) 


A 

1 

T 

L 

E 

S 

3 

2 

7 

Col umn 

i 

2 

3 

A 

5 

6 

7 

8 

9 

10 

1 1 

Field  1 

Field 

2 

Field  3 
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In  this  example,  variable  NAME  would  have  the  computer  representation 
of  ' T I TL 1 stored  in  it,  N would  contain  1 ES  "fet1  represents  a blank) 
and  J would  contain  the  alphanumeric  representation  of  '327'.  It  is 
important  to  note  that  this  "value"  of  J does  not  correspond  to  an  integer 
value  of  327. 

Mixed  Formats 


In  many  instances  integer,  floating  point,  and  alphanumeric  data 
all  appear  on  the  same  card.  It  is  simply  a matter  of  combining  the 
different  formats  in  the  correct  sequence  to  read  this  data.  Also, 
many  times  certain  columns  of  the  data  card  will  want  to  be  skipped. 
To  accomplish  this,  the  nX  Format  is  used  where  n is  the  number  of 
columns  to  be  skipped. 

Example : 

READ  103,  X,  I,  N , Y 
103  FORMAT  (F6.2,  AX,  12 , A3,  2X,  FA.1) 


Co  1 umns 


The  value  of  X will  be  A 7 • 3 , the  value  of  X wi 1 1 be  7 . the  "va 1 ue" 
of  N will  be  ABC,  and  the  value  of  Y will  be  672.2.  Notice  that  the  7 


i 
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which  appears  in  column  17  is  not  assigned  to  any  variable  since  it  is 
in  a column  which  is  skipped. 


For  each  item  of  data  that  is  needed  by  the 
READ  statement  and  FORMAT  statement  is  included, 
card  in  the  "User's  Manual." 


program,  the  associated 
along  with  a sample  data 


I 

I 


MAIN  PROGRAM  USERS  GUIDE 


1.  Required  Main  Data  and  Format 

The  only  data  required  by  MAIN  are  MASTER,  MFLAG,  and  EPSLON. 

If  NFLAG  is  "1"  the  program  stops  after  the  TOP  or  CORRCT  subroutine 
calculations  are  complete.  The  TOP  subroutine  calculates  the  flow,  mass, 

i 

and  pollutant  parameter  for  each  branch  of  the  system.  If  NFLAG  is  "0" 
the  program  goes  from  TOP  into  the  RM  subroutine  where  the  survey  is  planned. 

MASTER  is  the  largest  parameter  number  to  be  used  in  the  current  study. 
The  master  list  of  parameters  is  given  in  Table  II.  If  in  a given  study  no 
parameter  listed  above  total  phosphates  is  to  be  included  then  MASTER  is  18. 

If  EPSLON  = 0 the  CORRCT  subroutine  is  not  called.  For  EPSLON  = I 
the  CORRCT  subroutine  is  called. 

MASTER  and  MFLAG  are  read  in  with  the  following  statements: 

READ  (1,900)  MASTER 
READ  (1 ,900)  NFLAG 
READ  (1,900)  EPSLON 
900  FORMAT  (15) 

The  following  card  layout  is  for  MASTER=l8 

12  3 4 5 

The  following  car  layout  is  for  NFLAG=  1 or  for  EPSLON  = I 

12  3 4 5 


i i i i ii 
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II.  Changing  System  Planning  Program  Dimensions 

The  program  and  subroutines  have  been  written  so  that  only  cards  in 
the  MAIN  program  need  to  be  changed  to  redimension  the  system.  The  quan- 
tities that  can  be  varied  are:  the  number  of  pollutant  parameters,  the 

number  of  sources  (processes),  the  number  of  branches,  the  number  of 
pieces  of  equipment. 

The  specific  items  in  the  MAIN  program  which  must  be  changed  to  re- 
dimension the  system  are  given  in  the  following  paragraphs: 

To  handle  more  than  25  sources  (where  i = number  of  sources)  in  the 
MAIN  program, 

change  MSORS  = 25  to  MSORS  = i and  in  the  Dimension  statement 
the  arrays  to  modify  are: 


XNAME  (i) 
YNAME  (i) 

FLOW  (i) 

CHANGE  SPLIST  (i  ,25) 

CAP  (i) 

A (i ,25) 

TO  C ( i ,25) 

P (i ,25) 

SUM  (i) 

SUMA  (i) 


To  handle  more  than  25  branches  (where  n = number  of  branches)  in 
the  program, 

change  MBRNC  = 25  to  MBRNC  = n 


XNAME  (25) 
YNAME  (25) 

FLOW  (25) 
SPLIST  (25,25) 
CAP  (25) 

A (25,25) 

C (25,25) 

P (25,25) 

SUM  (25) 

SUMA  (25) 


and  in  the  DIMENSION  statement  the  arrays  to  modify  are: 


I 


A (25,25) 

A (25, n) 

Y (25,25) 

Y (n ,25) 

X (25,25) 

X (n ,25) 

EFF  (25) 

EFF  (n ,25) 

NBRNCH  (25,2) 

NBRNCH  (n  ,2) 

IBN  (51,25) 

IBN  (51,n+l) 

PC  (25,25,25) 

CHANGE 

PC  (25,25,n) 

AMAR  (25,27) 

AMAR  (25,  n+2) 

NALOW  (25,25) 

TO 

NALOW  (25, n) 

BRN  (25) 

BRN  (n) 

BRANCH  (25,25) 

BRANCH  (n , 25) 

NROUR  (25) 

NROUT  (n) 

NFLOW  (25) 

NFLOW  (n) 

To  handle  more  than  25  parameters  (m  = number  of  parameters) , 


the  MAIN  program, 


change  MPARM  = 25  to  MPARM  = m 


and  in  the  DIMENSION  statement  the  arrays  to  modify  are: 


SPLIST  (25,25) 

SPLIST  (25, m) 

C (25,25) 

C (25, m) 

Y (25,25) 

Y (25, m) 

X (25,25) 

X (25, m) 

P (25,25) 

P (25, m) 

NPLIST  (25) 

NPLIST  (m) 

POLN  (25,5) 

POLN  (m,5) 

NTEMP  (25) 

NTEMP  (m) 

SAMFRE  (25) 

CHANGE 

SAMFRE  (m) 

NMA  (25) 

NMA  (m) 

MENAME  (25,3,5) 

MENAME  (m,3,5) 

PCRM  (25,3) 

TO 

PCRM  (m,3) 

PMDATA  (25,3,l»,5) 

PMDATA  (m, 3,^,5) 

IDO  (25,25) 

IDO  (m,25) 

PM  (25,25,25) 

PM  (m,25,3) 

AMAR  (25,27) 

AMAR  (m,27) 

PC  (25,25,25) 

PC  (m,25,25) 

NALOW  (25,25) 

NALOW  (m,25) 

EQUSED  (25,75) 

EQUSED  (m,75) 

NSET  (25,3) 

NSET  (m, 3) 

FLGPT  (25) 

FLGPT  (m) 

RANK  (25) 

RANK  (m) 

BRANCH  (25,25) 

BRANCH  (25, m) 

ELE  (25) 

ELE  (m) 

To  handle  more  than  75  pieces  of  equipment  (j  = number  of  equipment) 


in  the  MAIN  program, 


change 

MEQ  - 75 

to 

MEQ  = j 

and  in  the 

DIMENSION  statement: 

change 

VNSP  (75) 

to 

VNSP  (j) 

SMEQTI  (75) 
EQUSED  (25,75) 
EU  (75) 

EQTIME  (75) 
EQNAME  (75,5) 

CHANGE 

TO 

SMEQTI  (j) 
EQUSED  (25, j 
EU  (j) 

EQTIME  (j) 
EQNAME  (j ,5) 

It  should  be  noted  that  these  are  the  only  changes  required  on  the 
array  dimensions.  There  is  no  need  to  make  modifications  in  the  sub- 
routines. Also,  it  is  very  important  that  the  array  parameters  MSORS, 
MBRNC,  MPARM,  MEQ  MBRP1,  and  MP2  and  the  array  subscripts  correspond 
exactly.  For  example,  if  MPARM  = 31,  then  every  array  that  is  affected 
by  an  increase  in  parameter  capability  must  have  "31"  appearing  in  the 
proper  subscript  position  as  explained  in  that  section.  Care  should 
also  be  exercised  that  every  array  has  its  proper  subscript  changed. 

Any  omission  may  lead  to  errors  that  may  not  be  immediately  apparent. 
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PROCESS  MODEL  USER'S  GUIDE 
. PROCESS  MODEL  INTRODUCTION 

The  process  model  is  used  to  define  a manufacturing  plant  in  terms 
of  the  pollutants  produced  by  the  manufacturing  processes  as  a function 
of  the  production  levels  of  the  various  plant  components.  In  order  that 
the  model  be  applicable  to  more  than  one  plant,  manufacturing  operations 
are  subdivided  on  the  basis  of  "objectives"  and  then  further  subdivided 
into  "processes"  required  by  these  objective.  "Objectives"  are  identified 
in  terms  of  general  plant  products  or  objectives.  A given  plant  may  have 
one  or  several  "objectives".  Examples  are:  continuous  TNT;  TNT  (batch); 

paint  stripping  of  artillery  shells.  The  processes  are  manufacturing 
operations  which  are  definable  in  terms  of  production  level  and  resultant 
pollutants.  The  processes  are  cperations  which  may  be  components  of  one 
or  several  "objectives"  and  may  occur  in  one  or  several  plants.  Models 
are  constructed  for  each  process  which  may  occur  in  any  of  the  plants  to 
be  surveyed.  A plant  process  model  is  obtained  by  assembling  the  appro- 
priate process  models.  The  input  to  this  model  is  the  production  level 
of  each  of  the  processes  making  up  the  plant  model.  The  output  is  an 
ordered  listing  of  the  types  and  levels  of  pollutants  generated  by  pro- 
ducing products  at  a certain  level.  Also  listed  are  the  points  in  the 
effluent  drain  system  at  which  the  pollutants  enter.  Any  time  varia- 
tion in  pollutant  flow  requires  a new  listing,  a listing  for  each  state. 

The  pollution  potential  of  an  effluent  is  measured  in  a number  of 
ways.  One  type  of  measurement  is  to  determine  the  concentration  of  a 
particular  element  or  compound.  An  example  of  this  type  of  measurement 
would  be  the  measure  of  sulfates  in  milligrams  per  liter.  Another  type 
of  measurement  measures  a characteristic  of  the  system  not  necessarily 
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0001 

0002 

0003 

0004 
0003 
0006 

0007 

0008 
0009 
0010 
001  1 
0012 

0013 

0014 

0015 

0016 

0017 

0018 
0019 
0020 
0021 


f 

I 

i 

l 


STNT 


SUBROUTINE  STNT  (FLOW,  SPLIST, 
DIMENSION  FLOW  (MSORS),  SPLIST 


FLOW  (JM)  - 

= 0.01 

* CAP 

SPLI 

1ST 

(JM 

,1)  = 

1.5 

SPL 

1ST 

(JM 

,3)  = 

40000 

SPLI 

1ST 

(JM 

,4)  = 

40 

SPLI 

1ST 

(JM 

,7)  = 

1000 

SPLI 

1ST 

(JM 

,8)  = 

8 

SPLI 

1ST 

(JM 

,9)  = 

20000 

SPLI 

1ST 

(JM 

,10)  = 

10 

SPLI 

1ST 

(JM 

, 12) 

= 2500 

SPLI 

1ST 

(JM 

,15)  = 

1500 

SPLI 

1ST 

(JM 

,16)  = 

5.0 

SPLI 

1ST 

(JM 

,18)  = 

0.5 

SPLI 

1ST 

(JM 

,19)  = 

10.0 

SPLI 

ST 

(JM 

, 20) 

= 20 

SPLI 

1ST 

(JM 

,22)  = 

0.0 

SPLI 

ST 

(JM 

,23)  = 

5.0 

SPLI 

ST 

(JM 

,24)  = 

5000 

RETURN 

END 


CAP,  JM,  MSORS, 
(MSORS,  MPARM) , 
0.10 


MPARM) 

CAP  (MSORS) 


EXAMPLE  PROCESS 


Figure  I 


17 


I 

I 

I 

t 

I 

I 

I 

I 

! 

I 

i 


related  to  only  one  chemical  element  or  compound.  Examples  of  this  type 
of  measurement  are  conductivity,  color,  and  pH.  To  speak  of  in  general 
of  the  numbers  that  define  the  level  of  pollution  potential,  whether  they 
are  of  one  or  the  other  of  the  above  types,  the  term  "pollution  para- 
meter" will  be  used . 

Figure  I i s an  example  of  a computer  program  subroutine  model  of 
a process.  Statement  0001  gives  subroutine  name.  Statement  0002  is  the 
DIMENSION  statement.  Statement  0003  computes  the  effluent  flow  from  this 
process  in  million  gallons  per  day  (MGD)  for  a given  operating  level, 

CAP  ( J M ) in  % capacity  for  JM.  The  remaining  lines  0004  through  0019 
give  the  pollutant  parameters  necessary  to  define  this  process.  These 
parameters  must  be  identified  by  a number  assigned  in  a master  list  given 
in  Table  II.  For  example,  SPLIST  (JM,3)  = 40000  gives  the  conductivity  of 
the  effluent  of  the  process  numbered  JM.  Three  is  the  parameter  number  for 
conductivity  from  the  master  list. 

Further,  each  process  is  identified  by  a mnemonic  from  a master  list. 

A list  is  given  in  Table  I. 

Finally,  the  process  modeling  is  complete  for  a given  plant  when  all 
of  the  processes  and  their  operating  levels  are  identified  for  that  plant 
and  the  computer  program  calls  in  the  appropriate  subroutines  which  defines 
these  processes.  This  information  includes  a listing  of  the  points  in  the 
drain  system  at  which  the  pollutants  enter. 
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II.  GENERAL  DESCRIPTION  OF  PROCESS  MODEL  PROGRAM 


The  process  model  is  identified  by  a mnemonic  name  for  a specific 
operation  within  a plant  and  is  defined  in  terms  of  production  level 
capacity  and  the  parameters  which  measure  the  pollutants  found  in  the 
effluent  of  that  operation.  A master  list  of  process  models  is 
available  in  Table  I.  These  models  are  contained  in  the  main  program 
card  deck. 

During  the  preliminary  survey  all  processes  which  make  up  the 
plant  operation  are  identified  in  terms  of  the  proper  mnemonic  from 
the  master  list.  The  expected  operating  capacity  in  % of  full  load  is 
noted  for  each  process.  The  mnemonic  and  capacity  are  placed  on  a 
data  card  for  read  in  by  the  PROCES  subroutine.  These  process  cards 
are  read  in  the  order  determined  by  the  topological  numbering  dis- 
cussed on  page  39  under  ordering  of  data  and  coding.  There  are  no 
restrictions  on  the  numbering  or  ordering  of  the  processes  except  that 
once  the  numbering  has  been  assigned  it  must  be  used  here  and  in  the 
topology  part  of  the  program. 

As  the  subroutine  PROCES  reads  in  the  process  data  cards  it  forms 
a plant  model  in  terms  of  a SPLIST  matrix.  The  SPLIST  matrix  has  as 
its  rows  the  different  processes  present  and  as  its  columns  the  parameter 
values,  the  parameters  in  order  as  defined  by  the  master  list  of  Table  II 
and  processes  in  the  order  read  in.  An  abbreviated  SPLIST  matrix  format 
example  for  the  first  five  process  models  as  well  as  the  first  four 
pollutant  parameters  is  described  for  BAAP. 
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ORDER  IN  WHICH  FOUND  IN  PROGRAM  AND  MNEMONICS 

PAPC  - AOP  Compression 

PAOP  - AOP  - Flow  3 

PSAC  - Sulfuric  Acid  Concentrator 

PNAC  - Nitric  Acid  Concentrator 

PBDN  - Blowdown  at  BAAP 

PBIX  - Ion  Exchange  at  BAAP 

PCOW  - Cooling  Water  - Flow  10 

PNBP  - Nitration  and  Boiling  Operations 

PBPP  - Beater  and  Poaching  Operation 

PBEX  - Benzene  Extraction 

PNCS  - Nitrocellulose  into  Ball  Power 

PNGS  - Ball  Power  with  Nitroglycerin 

PPND  - Pond  at  Head  of  Sewer 

PAPP  - AOP  - Flow  A 

PSWG  - Sewer  Influent  - Flow  into  STP  Process 

PCWZ  - AAP  - Cooling  Water  - Flow  11 

PPAS  - Purification  Acid  Scrubber 

PPCU  - Purification  Clean-Up 

PPCW  - Purification  Cooling  Water 

PFSW  - Finishing  Scrubber  Water 

PROCESSES 

(SOURCES) 

TABLE  1 


CORRESPONDING 

SUBROUTINES 

SAPC 
SAOP 
SSAC 
SNAC 
SBDN 
SB  I X 
SCOW 
SNBP 
SBPP 
SBEX 
SNCS 
SNGS 
SPND 
SAPP 
SSWG 
SCWZ 
SPAS 
SPCU 
SPCW 
SFSW 
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Table  1 Continued: 


ORDER  IN  WHICH  FOUND  IN  PROGRAM  AND  MNEMONICS  CORRESPONDING 

SUBROUTINES 

PFCU  - Finishing  Clean-Up  SFCU 

PFCW  - Finishing  Cooling  Water  SFCW 

PTCW  - Tank  Car  Wash  STCW 

PSAR  - Sulfuric  Acid  Recovery  Cooling  Water  SSAR 

PSAW  - Sulfuric  Acid  Recovery  Waste  Acid  SSAW 

PANP  - Acid  Neutralization  Plant  (Radford)  SANP 

PC  IX  - Ion  Exchange  at  Newport  - C 1 i ne  SC  I X 

PZBM  - Ion  Exchange  at  Newport  - A line  SZBM 

PCWB  - Power  House  Cooling  Water  at  Newport  SCWB 

PAPN  - AOP  at  Newport  SAPN 

PNDN  - NAC/DN  at  Newport  SNDN 

PGRB  - Grab  Sample  at  M-6  at  Newport  SGRB 

PSAZ  - Sulfuric  Acid  Recovery  at  Newport  SSAZ 

PTNT  - TNT  Lab  at  Newport  STNT 

PPHZ  - Pumphouse  tt 2 at  Newport  SPHZ 

PSTP  - Sewage  Treatment  Plant  Influent  at  Newport  SSTP 


PROCESSES 

(SOURCES) 


TABLE  I 
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Table 

PCWL  - 
PODW  - 
PCOA  - 
PCOB  - 
PRWS  - 
PAPS  - 
PADR  - 
PZPP  - 
PCAP  - 
PAPK  - 
PAHC  - 
PZPD  - 
PCAD  - 


1 Continued 


Cooling  Water  Blowdown 

SCWL 

Oi 1 /Descalar  Water 

SODW 

Cutter  Oil  - Line  A 

SCOA 

Cutter  Oil  - Line  B 

SCOB 

Rinse  Water  Supply 

SRWS 

Alkaline  Paint  Stripper 

SAPS 

A1 ka 1 i ne  De-Rus t 

SADR 

Zinc  Phosphate  Dip  - 

Paint  Stripping 

SZPP 

Chromic  Acid  Rinse  - 

Paint  Stripping 

SCAP 

Alkaline  Pre-Cleaner 

SAPK 

Alkaline  Hot  Cleaner 

SAHC 

Zinc  Phosphate  Dip  - 

Detrex 

SZPD 

Chromic  Acid  Rinse  - 

Detrex 

SCAD 

PROCESSES 


TABLE  1 
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1 . 

pH 

2. 

COD 

3. 

CONDUCTIVITY 

4. 

TOC 

5. 

AMMONIA  (NH^)  - N 

6. 

DISSOLVED  SOLIDS 

7. 

NITRITE/NITRATE  (N) 

8. 

TOTAL  N (KJELDAHL) 

9. 

TOTAL  SOLIDS 

10. 

SUSPENDED  SOLIDS 

1 1 . 

TOTAL  HARDNESS 

12. 

SULFATES 

13. 

MBAS  . 

14. 

SULFITES 

15. 

CHLORIDES 

16. 

VOL.  SUSP.  SOLIDS 

17. 

SODIUM 

18. 

TOTAL  PHOSPHATES 

19. 

TNT 

20. 

CALCIUM 

21 . 

COLOR 

22. 

ALKALINITY 

23. 

TURBIDITY 

24. 

ACIDITY 

25. 

MASTER 

BOD 

LIST  OF  PARAMETERS 

TABLE  II 
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PH 

COD 

CONDUCTIVITY 

TOC 

V-  Pol  1 utant  Names 

PPND 

6.0 

20.0 

0.0 

1 1.0 

in  Columns 

PAPC 

6.2 

20.0 

0.0 

1 1.0 

4 — i 

PAOP 

6.2 

20.0 

0.0 

1 1.0 

Pol  1 utant 
Parameters 
k — 1 

PAPP 

6.2 

20.0 

0.0 

11.0 

PNAC 

2.5 

25.0 

0.0 

22.0 

Process  Model  names  in  rows 

ABBREVIATED  SPLIST  MATRIX 
TABLE  I I I 


The  PROCES  subroutine  is  written  such  that  if  a process  mnemonic 
and  capacity  is  read  i n on  a data  card  and  there  is  no  library  sub- 
routine defining  this  process,  a row  will  be  reserved  in  SPLIST  with 
parameter  data  to  be  supplied  later. 

The  PROCES  subroutine  also  accepts  data  which  will  alter  one  or 
more  parameter  for  any  process  present.  If  the  user  feels  that  at  a 


I 


2k 


particular  plant  site  a different  parameter  value  gives  a better 
model,  the  parameter  can  be  changed.  Further  in  a more  extreme  case, 
all  new  parameters  can  be  read  in  for  a given  process  if  existing 
data  and  the  judgment  of  the  user  warrant  the  change.  Clearly,  this 
latter  capability  can  be  used  to  include  a process  in  a survey  planning 
when  the  process  is  not  in  the  process  library. 
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111.  GENERAL  DESCRIPTION  OF  PROCESS  MODELS 

| 

Flow  rates  in  MGD  are  based  on  the  production  capacity 
of  the  process  model  at  the  time  it  was  evaluated.  A summary 
of  each  process  model  as  of  July  1,  197**,  is  enclosed  by  loca- 
tion, code  name,  flow  rate,  capacity  at  1001,  capacity  at 
evaluation,  and  other  characteristics  to  clarify  the  process 

9 

ft 

description  as  deemed  necessary. 
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i 
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BAAP  Badger 
PAPC  (.57  MGD) 

PAOP  (.**8  MGD) 
PSAC  (1.6  MGD) 

PNAC  (.30  MGD) 

PBDN  (.08  MGD) 

PBIX  (.01  MGD) 

PC0W  (.30  MGD) 

PNBP  (.50  MGD) 
PBPP  (.20  MGD) 

PBEX  (.A5  MGD) 

PNCS  (2.0  MGD) 

PNGS  (.09  MGD) 
PPND  (.A0  MGD) 

PAPP  (.89  MGD) 


6 units  at  50  T/D.  Total  300  T/D.  Evaluated 
at  75  T/D.  3~ 8- 5 . A0P  Compression 

Same  as  PAPC  - (Cooling  water) 

3 units  at  150  T/D  each.  3“ 8- 5 • Evaluated 

at  110  T/D  (Cooling  water)  Sulfuric  Acid  Concentrator 

16  towers  at  30  T/D  each.  3-8-5.  Evaluated 
at  50  T/D  (Cooling  water)  Nitric  Acid  Concentrator 

A boilers  at  175,000  lb  stream/hr  each. 

Continuous  at  5 per  cent  of  production.  Blow  down 

7 columns  at  150,000  lb  stream/hr  each. 

Evaluated  at  A columns.  A00  GPM  for  10  min. 
backwash.  132  GPM  for  15  min.  slow  rinse. 

26A  GPM  for  50  min.  fast  rinse.  Once  every 
30  hours.  Ion  Exchange 

Cooling  water  from  boiler  at  25  percent 

capacity.  Continuous.  1.5x  stream  production.  Cooling  Water 

17.000  lb/day.  8 am  - 3 pm  weekdays.  Nitration  and  Boiling 

120.000  lb/day.  Evaluated  at  11,000  lb/day 
from  E-line.  Continuous.  Beater  and  Poaching 

610.000  lb/mo.  of  reprocessed  nitrocellulose. 

Intermittent.  Benzene  Extraction 

2.700.000  lb/mo.  of  ball  powder.  Evaluated 

at  700,000  lb/mo.  3-8-5  Nitrocellulose  into  Ball  Powder 

Same  as  PNCS. 

Pond  overflow  based  on  the  water  treatment  level. 

A5  MGD.  Evaluated  at  12  MGD.  Continuous.  Pond 

Same  as  PAPC  - (Process  water). 
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PSWG  (.27  MGD)  - Based  on  BAAP  facilities.  Continuous.  Sanitary  Sewer 

PCWZ  (.20  MGD)  - Cooling  water  from  compressors,  fans,  and 
pumps.  50  Percent  of  capacity  on  weekdays 
and  25  percent  on  weekends.  Continuous. 

NOTE:  All  values  for  BAAP  are  based  on  averages  in  the  assay 

or  data  sheets.  Additional  element  values  were  added 
to  the  assay  values  for  cooling  water  and  process  water 
effluents  by  using  the  data  sheets  for  raw  water  and 
process  water. 

RAAP  (The  £en  processes  for  one  50  T/D  continuous  TNT  line).  (Radford) 

PPAS  (.0015  MGD)  Purification  Acid  Scrubber 
PPCU  (.02  MGD)  Purification  Clean  Up 
PPCW  (.90  MGD)  Purification  Cooling  Water 
PFSW  (.0024  MGD)  Finishing  Scrubber  Water 
PFCU  (.01  MGD)  Finishing  Clean  Up 
PFCW  (.15  MGD)  Finishing  Cooling  Water 
PTCW  (.015  MGD)  Tank  Car  wash 

PSAR  (1.5  MGD)  Sulfuric  Acid  Recovery  Cooling  Water 

PSAW  (.01  MGD)  Sulfuric  Acid  Recovery  Waste  Acid 

PANP  (.135  MGD)  - unique  to  RAAP.  Acid  Neutral izat ion  Plant 
NAAP  Newport 

PCIX  (.3214  MGD)  - 3 cation  exchanges  with  2 operations  during 

evaluation.  6,000  gal.  tank  for  bulk  storage. 

Evaluated  at  90  percent  capacity.  Ion  exchange  Newport-A 

PZBM  (.0381  MGD)  - 2 anion  exchanges  with  1 operating  during 

evaluation.  Evaluated  at  90  percent  capacity. 

Ion  exchange  Newport-C 

PCWB  (.17  MGD)  - 120  GPM  normal  blowdown.  Evaluated  at 

90  percent  capacity.  Power  House  Cooling  Water  Newport 

PAPN  (.30  MGD)  - 350  T/D  of  65  percent  HN03.  100  percent  acid 

basis.  A0P  at  Newport 

- 2 NAC  units  with  100  T/D  capacity  of  98.5 

percent  acid.  4 denitiating  units  with  290  T/D 
of  TNT  waste  acid  each.  Nitric  Acid  Concent rator-Newport 


PNDN  (.30  MGD) 
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PGRB  (.002  MGD)  - Unique  to  NAAP  Grap  Sample  at  M-6  Newport 

PSAZ  (.30  MGD)  - Regeneration  of  7.4  T/D  of  40  percent  Oleum 

and  326  T/D  of  20  percent  Oleum.  Sulfuric  Acid  Recovery-Newport 

PTNT  (.10  MGD)  - 5 lines  at  50  T/D  each.  Evaluated  at  one  line 

with  outfall  unique  to  NAAP.  TNT  Lab-Newport 


PPHZ  (.02  MGD)  - Evaluated  at  90  percent  capacity.  Pumphouse  #2-Newport 


PSTP  (.50  MGD)  ■ Based  on  NAAP  facilities.  Sewage  Treatment  Plant  I nf 1 uent -Newport 


NOTE:  Values  were  based  on  conference  at  Clemson  University 

in  January  1974.  All  operations  are  continuous  and 
100  percent  capacity  was  used  unless  noted  otherwise. 


LAD  Letterkenney 

Shell  Production  - 160,000/mo.  Evaluated  at  100,000/mo. 

PCWL  (.0286  MGD)  * Sample  Point  J Cooling  Water  Blowdown 

P0DW  (.154  MGD)  * S-  nple  Point  K Oil/Descaler  Water 

PC0A  (.01  MGD)  - G-12  (10,000  gal.  cap.  tank)  Cutter  Oil  Line  A 

PC0B  (.01  MGD)  - G-I3  (10,000  gai.  cap.  tank)  Cutter  Oil  Line  8 
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Paint  Stripping  - 3 to  5 percent  of  the  shells  are  rejected. 

Capacity  is  500-700  shells/day/shift.  There 
are  2 lines  of  which  1 line  was  used  during 
evaluation  not  at  fall  capacity.  Continuous 
overflow  for  rinse  is  2 GPM. 

PRWS  (.0056  MGD)  Rinse  Water  Supply 

PAPS  (.00035  MGD)  Alkaline  Paint  Stripper 

PADR  (.00035  MGD)  Alkaline  Derust 

PZPP  (.00035  MGD)  Zinc  Phosphate  Dip 

PCAP  (.00035  MGD)  Chromic  Acid  Dip 

Point  Preparation  - "DEXTREX"  process  of  batch  dumps  which 

last  10  minutes  at  the  end  of  each  work 
period.  Continuous  overflow  for  pre- 
cleaner and  hot  cleaner  A GPM. 

PAPK  (.001  MGD  Weekly)  Alkaline  Pre-Cleaner 


PAHC  (.0012  MGD  Monthly)  Alkaline  Hot  Cleaner 

PZPD  (.0011  MGD  Monthly)  Zinc  Phosphate  Dip-  Detrex 

PCAD  (.0009  MGD  Daily)  Chromic  Acid  Dip-Detrex 

Notes : Averages  of  days  2,  6,  and  9 were  used  for  point  "J". 

Averages  for  days  8-13  were  used  for  point  "K". 


I V . REQUIRED  PROCESS  SUBROUTINE  DATA  INPUT  AND  FORMAT 


The  following  paragraphs  describe  the  data  required  by  the 
PROCES  subroutine  and  the  card  format  required  by  the  subroutine. 

The  cards  and  data  will  be  discussed  in  the  order  in  which  they  must 
be  read  in.  The  data  required  is: 


1.  The  number  of  processes  in  the  plant  model,  N 

2.  The  process  names,  XNAME  (JM) , and  the  percent 
of  full  production  capacity,  CAP  (JM)* 

3.  A flag  indicating  whether  or  not  there  are  to 
be  pollutant  parameter  modification  for  any  of 
the  processes,  NOD 

A.  If  there  are  modifications,  give  the  number 
of  parameters  to  be  changed,  NMOD 

5.  For  each  parameter  to  be  changed,  give  the 
process  name,  the  pollutant  number,  and  the 
new  pollutant  value.  YNAME  (NM ) , ELE  (NM)  , 

VAL  (NM)  - One  card  for  each  parameter  to  be 
changed . 

6.  A flag  indicating  whether  or  not  there  are  to 
be  any  complete  processes  changed,  NPRO 
(Could  be  a new  process) 

7.  If  there  are  processes  to  be  changed,  give 
the  number  of  processes  to  be  changed,  NUMP 

8.  Read  in  the  number  of  parameters  currently  in 
the  master  list,  NELE . 

9.  For  the  first  process  to  be  changed,  read  in 
the  process  name  and  the  % of  full  capacity 
flow.  ZNAME  (JN)  , FLW  (JN) 

10.  One  card  for  each  parameter  In  the  master  list 
in  the  same  order  as  in  the  master  list,  giving 
the  new  parameter  values,  SPLIST  (JE,  JK) 

11.  Repeat  9 and  10  for  each  process  to  be  modified. 


*Must  be  readin  in  the  same  order  as  topological  numbering. 


30 


The  data  card  formats  are  read  in  statements  for  the  PROCESS 
data  and  are  now  given. 


1.  The  number  of  processes  in  the  plant  model  N.  Example 


for  16  proceses: 


READ  (1 , 1 0) N 
10  FORMAT  (IX,  12) 


Card  Columns 


2.  One  card  for  each  process  giving  process  mnemonic  and  the  % 

of  full  production  capacity.  Example:  PTNT  at  500%  capacity. 

READ  (1,200)  XNAME  (JM)  , CAP  (JM) 

200  FORMAT  (IX,  A4 , IX,  F6.2) 

1 2 3 4 5 6 7 8 9 10  1 1 12  <c—  Card  Columns 
f|  I P I T I N It  I I 5 I 0 |o  I . I 0 I 0 I 


3.  Flag  to  indicate  pollutant  parameter  modifications.  If  there 
are  modifications  N0D=1.  If  there  are  no  modifications  N0D=0. 
Example:  N0D=l 

READ  (1,400)  NOD 
400  FORMAT  ( 1 X , II) 


1 2 3 


Card  Columns 
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If  there  are  modifications  the  next  card  gives  the  number  of  parameters 
to  be  modified.  Example:  21  (bypassed  if  N0D=0) 

READ  (1,700)  NMOD 
700  FORMAT  (IX,  13 ) 

12  3^5  ^ Card  Columns 


5.  In  this  group,  one  card  for  each  parameter  to  be  changed.  Give 

the  process  mnemonic,  the  pollutant  number  and  the  pollutant  parameter 
value.  Example:  PPND,  pollutant  number  2,  pollutant  parameter  value 

21.15  (number  2 is  COD) 

READ  (1,750)  YNAME  (NM) , ELE  (NM) , VAL  (NM) 

750  FORMAT  (IX,  IX,  12,  IX,  F9.2) 


6.  This  is  a flag  card  to  indicate  if  any  complete  processes  are  to  be 

changed.  If  there  are  modifications  NPR0=1 . If  there  are  no  modifications 
NPR0=0 . Example:  NPR0=1 

READ  (1 ,875)  NPR0 
875  FORMAT  (IX,  ll) 

1 2 3^5  4 Card  Columns 
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7.  If  there  are  complete  processes  to  be  modified,  this  card  gives  the 

number  of  processes  to  be  changed.  Example:  12 

(bypassed  if  NPR0=0) 

READ  (1 ,890)  NUMP 
890  FORMAT  (IX,  12) 

Card  Columns 

8.  Next  read  in  the  number  of  parameters  currently  in  the  master  list. 

NELE  Examp1e:2A 

READ  (1 ,890) NELE 
890  FORMAT  (IX,  12) 

1 2 3 A 5 6 Card  Columns 


9-  For  the  first  process  to  be  changed,  read  in  the  process 
mnemonic  and  the  % full  load  capacity. 

Example:  ZNAME  (JN)  = PPND.FLW  (JN)  = 80 

READ  (1,895)  ZNAME  (JN)  , FLW  (JN) 

895  FORMAT  (IX,  AA , IX,  F6.A) 


1 2 3 A 5 6 7 8 9 10  11  12  13  1 A 15 

t 

| I P IP  |N  | D | 1 8 1 0 | . J 0 I 0 1 | | | | 


Card  Columns 


f 


I 
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10.  This  next  group  of  cards  (NELE  in  number),  one  for  each  parameter 
in  the  master  list,  gives  the  value  for  each  parameter  defining 
the  process.  These  must  be  in  the  same  order  as  in  the  master 
list  of  parameters . 

Example  card:  the  third  card  - conductivity  = 1800. 

DO  897  JK-1 , NELE 

READ  (1,898)  SPLIST  (JE.JK) 

898  FORMAT  (IX,  F9.2) 

897  CONTINUE 


123456789  10  Card  Columns 

/ 

I I 1 | 8 |o  | 0 | . | 0 | 0 | 


1 


11.  Cards  described  in  7 through  10  are  bypassed  if  NPRO  = 0. 

Procedures  described  in  9 and  10  are  repeated  for  each  process 
to  be  modified.  Clearly  this  procedure  can  be  used  to  include 
processes  in  a survey  plan  when  these  processes  are  not  included 
in  the  process  library.  However,  this  use  does  not  place  the 
new  process  in  the  library.  The  next  section  explains  how  a 
new  process  may  be  added  to  the  library. 


3^ 


* 


V.  Procedure  for  Adding  a Process  Model  to  the  Library: 

A.  Add  a DATA  statement  to  the  existing  set  of  data  statements 
to  initialize  the  new  process  name. 


DATA  P J_  2_  3AHP  1_  2 3/ 


B.  Add  an  IF  statement  to  the  existing  set  of  IF  statements  to 
allow  the  main  program  to  call  the  new  process 

I F (XNAME  ( JM) . EQ.  P J_  1_  3_  CALL  S J_  2 3 

New  Subroutine 

Process  Name  for 

Name  k Subroutine 

Characters  that  Defines 

the  New  Process 

Positions  1 2 3 being  the  same. 

C.  Add  a subroutine  with  the  following  form  to  define  the  new  process. 
SUBROUTINE  S _l_  2 3 (FLOW,  SPLIST,  CAP,  JM,  MSORS,  MPARM) 

DIMENSION  FLOW  (MSORS),  SPLIST  (MSORS .MPARM) , CAP  (MSORS) 

FLOW  ( JM)  = 0.01  * CAP(JM)  * (Some  value  to  represent  the  flow 
of  the  process  at  100%  capacity) 

SPLIST  ( J M , 1 ) = (value  of  first  element  in  process) 

SPLIST  (JM,2)  = (value  of  second  element  in  process) 


. These  are  values  of  pollutant  I through  N 

SPLIST  (JM,N)  = (value  of  N^h_  element  in  process1 

RETURN  The  SPLIST  array  must  be  redimensioned 

END  if  N > 25. 

If  a pollutant's  value  is  not  known  or  is  zero,  then  the  statement  de- 
fining that  pollutant  may  be  omitted  and  the  program  will  insert  a zero 
into  the  SPLIST  array  for  that  pollutant. 


35 


V | . Summary : 

The  user  must  be  aware  of  the  most  recent  master  lists 
of  process  models  and  pollutant  names  contained  in  the  system 
library  so  that  he  can  plan  his  overall  data  gathering  require- 
ments during  the  preliminary  survey.  He  must  not  feel  that  each 
new  operation  is  unique,  requiring  the  development  of  data  for 
the  addition  of  new  process  models. 

As  more  process  models  are  added  to  the  library,  the 
emphasis  in  preliminary  surveys  can  be  placed  on  the  determination 
of  the  flow  rate  at  full  production,  the  production  operating 
capacity  during  measurement  periods;  and  a few  measurements  of 
operations  similar  to  the  process  models  in  the  library  to 
verify  previously  recordec  data.  When  there  is  a large  de- 
viation in  a pollutant  parameter  the  methods  and  procedures 
of  the  operation  which  contributed  the  unfavorable  parameter 
should  be  investigated  and  compared  to  the  favorable  operation 
to  possibly  recommend  improvement  rather  than  merely  record 
the  unfavorable  parameter  as  a modification.  Confidence  in 
the  validity  of  the  data  will  increase  as  the  surveys  improve 
i n qua  1 i ty . 
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TOPOLOGY  MODEL  USER'S  GUIDE 
I.  Collection  of  Data: 

The  topology  data  of  a site  consists  of  the  following: 

1.  A map  of  the  sewer  lines  (include  flow  direction). 

2.  Identification  of  all  processes  (sources)  noting  where 
these  processes  empty  into  the  sewer  system. 

3-  Identification  of  all  treatment  facilities  and  where 
these  treatments  occur  in  the  sewer  system. 

*4.  Note  the  location  of  all  manholes,  weirs  or  other 
possible  sample  points  in  the  sewer  system. 

5.  Identification  of  pollutant  parameters  to  be  measured. 

6,  Settling  pond  areas  if  any-type  of  neutralization  if  any-type  of 
domestic  treatment  if  any 

figure  II  a and  II  b give  the  topology  data  for  Badger  AAP 
and  serves  as  a good  example  of  how  to  record  the  necessary  data. 

The  sewers  are  represented  as  directed  line  segments  showing  the 
flow  directions  and  interconnections. 

The  processes  are  represented  as  rectangles  and  are  identified  by 
number  and/or  some  meaningful  acronym.  The  location  of  the  point  where 
each  process  discharges  its  effluent  into  the  sewer  is  also  known. 

Treatment  facilities  are  represneted  by  ovals  and  are  identified 
by  an  acronym  describing  the  type  of  treatment.  (Acid  or  caustic  neu- 
tralization, settling,  domestic  sewage  treatment.) 

The  location  of  all  manholes,  weirs,  or  other  possible  sample 
points  are  represented  as  circles  on  the  sewer  line.  Further  iden- 
tifying and  describing  information  can  be  included  if  necessary  as  shown 
by  the  collumn  on  the  right  side  of  Figures  I a and  I b. 
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II.  Ordering  of  Data  and  Coding 

All  of  the  data  must  be  ordered  (numbered)  and  coded. 

All  processes  (sources)  must  be  identified  on  the  master 
list  for  which  subroutine  exist  in  the  porcess 
model.  This  master  list  is  given  in  Table  I.  Each  process 
is  then  assigned  the  proper  mnemonic  listed  in  the  left 
column  of  Table  I.  This  correct  coding  is  necessary  if  the 
program  is  to  call  the  correct  subroutine  (named  in  the  right 
hand  column)  which  defines  that  process.  Note:  In  the  sec- 
tion describing  "Processes,"  it  is  explained  how  a process 
in  the  master  list  can  be  modified  to  fit  the  needs  of  the 
user. 

After  the  processes  have  been  assigned  their  appropriate 
mnemonic,  the  processes  must  be  numbered.  The  numbering  order  is  ar- 
bitrary but  when  assigned  it  must  be  used  consistently  through- 
out the  program  (in  both  "Process"  and  "Topology"). 

Next  one  needs  to  identify  and  order  the  measure  points 
and  identify  the  flows  of  all  processes  which  pass  through 


each  measure  point.  This  ordering  and  identifying  is  best 
done  by  using  a linear  graph  as  shown  in  Figure  III  for  the 
Badger  AAP  system  of  Figure  II  a and  II  b. 


I 
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2.  Each  manhole  or  other  measure  point  is  represented  by 
a directed  line  segment  on  the  graph. 

3.  Each  leg  of  the  sewer  system  that  does  not  contain  a 
manhold  or  measure  point  is  contracted  to  "zero"  length. 

4.  The  direction  of  flow  in  each  line  segment  is  indicated 
by  an  arrow  head. 

5.  Partial  treatments  are  represented  by  ovals  containing 
an  identifying  code  and  a line  from  the  oval  to  the 
point  in  the  system  where  the  flow  is  treated. 

6.  The  essential  information  on  the  linear  graph  includes 
the  identity  and  ordering  of  all  processes,  identity 
and  ordering  of  measure  points,  and  flows  from  each 
process  which  pass  through  each  possible  measure  point. 

As  a rule  of  thumb,  one  can  begin  the  construction  of 

the  system  linear  graph  by  going  to  the  "preliminary  survey" 

data  such  as  shown  in  Figures  II  a and  II  b and  start  by  listing 

a process  approximately  the  greatest  number  of  measure  points 

away  from  the  outfall.  In  the  example  this  is  the  "overflow 

from  pond",  its  mnemonic  is  PPND.  We  shall  assign  this  as 

process  number  one  (1)  and  note  that  the  flow  enters  the 

system  at  a manhole  which  we  will  identify  as  measure  point 

number  one  (1)  and  represent  as  a directed  line  segment  on 

our  graph  in  Figure  IV. 


Beginning  Construction  of 
System  Linear  Graph 
(Step  n 
Figure  IV 
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Next  as  we  proceed  from  measure  point  number  one  to  the 
next  measure  point  to  which  we  assign  number  two,  note 
that  flow  from  the  process  identified  as  "AOP  Compression 
House"  must  be  included  at  measure  point  two.  The  process 
is  identified  by  its  mnemonic,  PAPC , and  assigned  the  number 
two,  (2).  In  Figure  IV  the  flows  from  processes  1 and  2 are 
included  at  measure  point  number  2 represented  by  a directed 
line  segment 

Figures  I a and  I b show  no  manholes  or  other  measure 
points  in  the  sewer  system  from  the  point  marked  NF  on 
Fig.  II  a to  the  point  marked  18  on  Fig.  II  b.  Thus,  all  these 
sections  of  sewer  line  are  contracted  to  a single  point  and 
the  Figure  V shows  the  linear  graph  including  measure  point  18. 


Construction  of  System 
Linear  Graph 
(Step  2) 

Figure  V 


It  is  noted  from  Figures  II  a and  II  b that  several  other 
flows  are  included  at  measure  point  18.  In  Figure  VI  several 
of  these  are  added  and  Figure  III  shows  the  completed  linear 
graph.  Figure  VI  shows  the  addition  of  processes  PAOP  (3), 
PAPP  (6),  PSAC  (A),  and  PNAC  (5).  Measure  points  A,  5,  6,  7, 
and  8 are  included.  Acid  neutralization  treatment  is  in- 
cluded as  ANP  oval  acting  on  flow  entering  measure  point  8. 


Construction  of  System 
Linear  Graph 
(Step  3) 

Figure  VI 
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There  are  restrictions  in  numbering  the  measure  points  (branches  of  the 


linear  graph).  The  restrictions  are: 


I.  The  numbering  must  start  at  1 and  continue  consecutively 


to  the  highest  number. 


2.  The  numbering  must  be  such  that  in  tracing  from  any  branch  j 


to  the  outfall  only  higher  numbered  branches  will  be  encountered. 


3.  Thus,  the  outfall  branch  must  have  the  highest  number. 


If  the  sewer  system  is  in  separate  parts  and  hence  has  multiple 


outfalls,  the  outfall  elements  of  the  separate  parts  are  connected  together 


and  a fictitious  new  outfall  element  is  included  in  the  linear  graph.  In 


the  input  data  a data  item  designated  NTOP  gives  the  program  the  information 


necessary  to  handle  this  situation.  Figure  VII  illustrates  the  addition 


of  the  fictitious  outfall. 


21  - outfall  Part  1 

22  - outfall  of  Part  2 

23  - outfall  of  Part  3 

24  - fictitious  outfall 


Multiple  Outfalls 
Figure  vil 


Construction  of  the  Topology  Matrix 


A convenient  way  of  tabulating  the  flow  information  for  entering  into 


the  computer  program  is  to  put  the  information  in  a matrix  form.  The  infor- 


mation required  for  each  branch  of  the  linear  graph  (each  possible  measure 


point)  includes: 


1.  Does  the  flow  of  source  i flow  through  branch  j? 

2.  Does  the  flow  undergo  any  treatment  in  branch  j?  If  so, 

what  treatment? 

A matrix  is  constructed  with  a column,  j,  for  each  branch  of  the 
system  and  a row,  i,  for  each  source  of  the  system.  The  entries  for  the 
i1^1  row  are  obtained  by  starting  at  the  if^  process  on  the  linear  qraph 
and  tracing  a path  through  the  branches  to  the  outfall.  If  the  path  does 
not  include  the  branch  the  ij  entry  is  zero.  If  the  path  does  include 

the  j ^ branch  an  appropriate  number  is  entered  in  the  ij  position.  That 
appropriate  number  depends  upon  whether  or  not  the  flow  undergoes  treat- 
ment in  that  branch  and  specifically  which  treatment.  If  the  flow  of 
source  i passes  through  branch  j and  no  treatment  occurs  at  branch  j, 
then  the  ij  entry  is  1.  If  the  flow  of  source  i passes  through  branch  j 

and  the  composite  flow  in  branch  j is  treated  at  branch  j,  then  the  ij 

entry  is  2,3,  or  A depending  on  the  type  of  treatment  as  follows: 
acid  or  caustic  neutralization  -2 

sett  1 i ng  pond  ~3 

domestic  sewage  treatemnt  -A 

Figure  VIII  gives  the  topology  matrix  for  the  system  of  Figure  III. 

IV.  Identification  of  Pollutant  Parameters 

It  is  necessary  to  supply  information  to  the  programs  which  tell 
the  particular  pollutant  parameters  that  are  to  be  measured  on  a survey. 
Each  of  the  processes  are  defined  in  terms  of  2 5 parameters.  However, 
all  of  these  may  not  be  measured  in  a given  survey.  These  25  parameters 
are  listed  and  numbered  in  a master  list  given  in  Table  II.  These  are  in 
the  same  order  used  in  the  process  model.  The  table  numbers  for  specific 
parameters  is  supplied  to  the  program  on  a data  card.  Also,  another  set 


of  cards  is  used  to  read  in  the  names  of  the  parameters  to  he  measured. 
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V.  Names  of  System  Branches 

The  program  was  written  to  require  names  as  well  as  numbers  to  be 
assigned  to  the  branches.  Up  to  an  eight  letter  mnemonic  is  used  here 
Meaningful  names  assigned  here  will  allow  more  rapid  interpretation  of 
print  outs. 


VI.  Treatment  Facilities  Data 

The  TOP  subroutine,  in  calculating  the  pollutant  valves  for  each 
branch  of  the  system,  takes  into  consideration  any  treatment  facilities 
in  the  system.  Treatment  types  modeled  include  acid  or  caustic  neutra- 
lization, settling  pond  (industrial),  and  domestic  sewage  treatment. 

Data  required  to  identify  and  define  the  types  of  treatment  used 
f o 1 1 ov/s  . 


1.  Clarifier  or  settling  pond  area  in  acres  (industrial  treatment) 

2.  Acid  or  caustic  neutralization  (industrial  treatment) 


MNEUT  / 1 
2 


<! 


ac i d neutra 1 i zat i on 
caustic  neutralization 


If  MNEUT  = 1 then  for  neutralizing  agent 


MENU 


Ca ( OH) 2 
CaC03 
Na2C03 
NaOH 


If  MNEUT  = 2 


NCAUS  f 1 

\ 2 


H2S04  , HCI  or  HN03 
C02 


If  NCAUS  = 1 


MACY 


1 H2S04 

2 HCI 

3 HNO3 
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3.  Domestic  sewage  treatment. 

ITYPE  = 0 trickling  filter 

= 1 activated  sludge 

TRICKLING  FILTER 

AREAP  primary  clarifier  area  in  acres 
TVOL  trickling  filter  volume  in  acre-feet 
TEMP  wastewater  temperature  (°C) 

R2  recycle  from  filter  effluent  (MGD) 

R3  recycle  from  secondary  clarifier  (MGD) 

Rk  recycle  to  filter  (MGD) 

TFAREA  trickling  filter  area  in  acres 
K20  = 0.23  for  1"  rock  media 

= 0.13  for  2V  rock  media 

AREAS  secondary  clarifier  area  in  acres 
ACTIVATED  SLUDGE 
AREAP  as  above 
AREAS  as  above 

ABVOL  aeration  basin  volume  in  (MG) 

MLSS  mixed  liquor  suspended  solids  (Mg/1) 

YN  net  yield  from  plant  log  data 

R1  activated  sludge  recycle  (MGD) 

VII.  Flagg  i ng 

The  flagging  information  required  by  the  program  is  generally  not 
data  collected  on  the  preliminary  survey  but  rather  is  data  based  on  the 
judgement  of  the  survey  officer.  The  judgement  may  be  based  on  experience 
but  will  also  be  based  on  results  of  trial  runs  of  the  planning  programs 
during  the  survey  planning  operation.  If  a flag  is  associated  with  a 


U3 


particular  pollutant  at  a particular  point  in  the  system,  measurement 
will  be  made  of  that  pollutant  at  that  point  and  resources  will  be 
committed  before  assignments  are  made  by  the  Resource  Matching  Model. 

Three  types  of  flagging  procedures  are  available. 

1.  Automatic  flagging. 

2.  Manual  flagging  of  sources 

3.  Manual  flagging  of  branches 

For  automatic  flagging  a subroutine  is  included  which  will  check 
the  total  mass  of  any  pollutant:  discharged  at  the  outfall.  The  program 
will  then  automatically  flag  any  source  which  contributes  more  than  X 
(decimal  fraction)  to  this  total.  The  source  is  flagged  for  the  measure- 
ment of  this  parameter  only.  The  survey  officer  can  specify  the  value 
of  X for  each  parameter. 

The  subroutine  also  permits  the  survey  officer  to  flag  for  measure- 
ment any  particular  source  for  any  particular  parameter  or  any  particular 
branch  for  any  particular  parameter. 

VIII.  Required  Topology  Subroutine  Data  Input  and  Format 

The  input  data  required  by  the  topology  part  of  the  survey  planning 
program  are: 

1 . The  number  of  sources , NS 

2.  The  number  of  branches,  NB 

e.  The  number  of  pollutant  parameters,  NP 

A.  The  number  of  outfalls,  NTOP 

5.  The  number  designation  from  Table  II  of  the  parameters  to  be 
measured  on  the  survey,  NPLIST  (l) 

6.  The  topology  matrix  data,  A (j.l) 

7.  The  branch  names,  NBRNCH  (l,J) 


1 

I 

f 
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*8.  The  length  of  the  survey  in  days  (LENGTH) 

*9.  The  number  of  noncompeting  parameters  ( I STOP) 

10.  The  parameter  names  POLN  (N,l) 

*11.  The  number  of  measure  methods  for  each  parameter  (NMA(N)) 

*12.  The  sampling  frequency  per  day  (SAMFRE(N)) 

13-  Sewage  treatment  data 

14.  The  flagging  levels  FlGPT(l) 

15.  The  source  and  parameter  manual  flags  NF,  NX(I),  NY ( I ) 

16.  The  branch  and  parameter  manual  flags  NFB,  NX  ( I ) , NY  ( I ) 

•'These  data  are  read  in  here  by  subroutine  CHKI  but  are  used  in 

Resource  Matching. 

The  data  cards  for  these  data  arc  reed  in  the  order  listed  above 
immediately  after  the  process  model  data  cards.  The  card  formats  are 
described  below:  (this  data  is  called  for  by  topology  subroutines.) 

The  NS,  NB , NP,  NTOP  are  read  in  on  one  card,  the  first  card,  with 
four  fields  each  five  long  and  right  justified. 


12345b78 


I I 1 • 61 


I 1 12  '1  | I 1 1 1 I o| 

f'  f 

\!  D s M D 


I I h 


NT01 


FORMAT  (4l5 


The  number  designations  from  Table  II  of  the  parameters  to  be  measured 
are  read  in  on  one  card,  with  25  fields  each  three  long  and  right  justi- 
fied. These  are  read  in  as  NPLIST.  The  following  card  layout  is  for 
measurement  of  parameters:  1,  3,  5.  6 etc. 


- ^ 


I 
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FORMAT  (2513) 

The  topology  matrix  is  read  in  next,  one  row  of  the  matrix  per  card 
and  the  rows  in  the  same  order  as  in  the  topology  matrix.  The  cards 
are  read  in  with  the  following  statements: 

DO  2 J i 1 , NS 

READ  (1,3)  ( A ( J , I ) , 1 = NB)) 

3 FORMAT  ((8011)) 

2 CONTINUE 

As  an  example,  the  following  card  layout  is  for  the  first  row  of  the 
topology  matrix  of  Figure  VI  II: 

1 2 3 4 5 6 7 8 9 10  11  12  13  14  15  16  17  18  19  20  21 

/ 110000000  0 0 0 0 0 0 0 0 1 2 0 3 


The  branch  names  are  read  in  an  order  with  the  following  statements: 
DO  30  I » 1 , NB 
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30  READ  (1,  901)  (NBRANCH  (I,J),  J = 1 ,2l 
901  FORMAT  (2A4) 

For  example,  branch  #1  may  be  assigned  the  name  P0ND1 . The  card 
format  would  be  as  shown  below.  One  card  is  used  for  each  branch, 
(up  to  8 characters  can  be  used) 

12  3 4 5 6 7 8 


The  length  of  the  survey  design  in  days  is  read  in  at  this  point, 
by  CHK  1 subroutine.  The  information  is  used  in  the  resource 
matching  subroutine.  The  sampling  frequency  and  constraints  are 
given  on  a per  day  basis.  If  survey  periods  different  from  a 
per  day  basis  need  to  be  considered  in  planning  this  entry  allows 
such  periods  to  be  considered. 

The  number  of  days  is  read  in  with  the  following  statements: 

READ  (1 , 903)  LENGTH 
903  FORMAT  (12) 

For  a one  day  period  the  card  layout  is  as  follows: 

1 2 3 4 5 6 7 8 


The  number  of  noncompeting  parameters  is  read  in  at 
this  point  by  CHKt.  The  information  is  used  in  the  resource 
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matching  subroutine.  The  program  is  written  so  that  parameters  which 
have  only  one  measurement  method  available  can  be  listed  first.  Avail- 
able resources  are  assigned  for  these  parameters  first.  If  an  ISTOP 
number  is  entered  to  designate  the  number  of  these  parameters,  then 
when  the  program  is  later  searching  for  alternatives  it  will  not 
search  back  through  these.  The  designation  is  ISTOP  and  it  is  read 
in  with  the  following  statements: 

READ  (1,100)  ISTOP 
100  FORMAT  (12) 

And  for  2 noncompeting  parameters  the  card  layout  follows: 


The  names  of  all  parameters  in  the  master  list  (Table  ll)  must 
be  listed  in  the  same  order.  The  program  saves  the  names  of  those 
parameters  previously  designated  to  be  included  in  this  survey.  One 
parameter  is  listed  per  card.  On  the  same  card,  for  each  parameter 
the  number  of  measure  methods  available  (NMA  (N) ) and  the  sampling 
frequency  (SAMFRE  (N))  for  that  parameter  are  listed.  The  NMA(N) 
and  the  SAMFRE(N)  are  used  in  the  resource  matching  subroutine  but 
are  read  in  here  by  subroutine  CHKI  while  entering  the  parameter 
name.  The  significance  of  the  NMA(N)  is  discussed  in  the  resource 
matching  section  of  the  manual.  The  SAMFRE(N)  is  the 


number  of  samples  for  day  per  measure  point  for  the  parameter. 

The  data  card  layout  is  according  to  the  format  statement: 

FORMAT  (5AA,  5X,  II,  Ax,  F6.1) 

and  for  the  parameter  alkalinity,  2 measure  methods,  and  A samples  per 
day  per  measure  point: 


123A56789IO 20 25  26  30 3JL 

1A1L1K1A1L1I1N1I  1T1Y1  1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 2 1 1 1 1 1 I 1 I 141.1 

POLN  (N) 

NMA(N)  SAM  FRF.  (N) 

The  sewage  treatment  data  is  read  in  next.  The  reading  format  and 
the  programs  using  the  data  are  written  so  that  all  cards  listed  must  be 
included.  Appropriate  data  is  included  on  the  cards.  When  a particular 
treatment  is  not  included,  the  cards  normally  carrying  this  data  are  left 
blank,  (i.e.  zeros  are  read  in)  This  data  is  read  in  by  the  TOP  subroutine. 

The  area  of  the  settling  pond  or  clarifier  is  read  in  (for  industrial 
treatment).  This  area  is  in  acres  and  is  read  in  with  the  statements: 

READ  20,  AREA 
20  FORMAT  (F5.2) 

The  card  layout  for  an  8 acre  pond  is 

1 a 3 h f <0  7 

/ 1 IS!  •toio\  I 1 


If  the  treatment  is  for  industrial  waste  and  involves  acid  or  caustic 
neutralization  the  next  card  supplies  the  value  for  MNEUT. 
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MNEUT  = 1 for  acid  neutralization  and 
MNEUT  = 2 for  caustic  neutralization 

MNEUT  is  read  in  with  the  following  statements 
READ  21 , MNEUT 
21  FORMAT  (II) 

The  card  layout  for  MNEUT  = 2 is 

/ixi — 

/ill  II  I 

If  the  treatment  is  acid  neutralization,  then  MENU  indicates  the 
neutralizing  agent. 


Ca (OH) 2 - 1 
CaC03  - 2 
Na  C03  - 3 
NaOH  - 4 


MENU  is  read  in  with  the  following  statements: 

READ  22,  MENU 
22  FORMAT  (II) 

The  card  layout  for  MENU  = 2 is 

/ ^ 3 v r 

Maui  I I 

If  the  treatment  is  caustic  neutralization,  then  NCAUS  indicates  the 
neut  ra I i z i ng  agent 

H2S04,  HCt  or  HNO3  - 1 
C02  - 2 
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NCAUS  is  read  in  with  the  following  statements: 

READ  23,  NCAUS 
23  FORMAT  (11 ) 

The  card  layout  for  NCAUS  = 1 is 

>73  Y 5 

/mTTT7 

If  NCAUS  = 1 then  MACY  gives  the  neutralization  agent 

H2S04  - 1 = MACY 

HCI  - 2 = MACy 

HN03  - 3 = MACY 

MACY  is  read  in  with  the  following  statements 
READ  2A,  MACY 
2 A FORMAT  (II) 

The  card  layout  for  MACY  = 3 is 


/ a 3 y r 
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If  the  treatment  is  for  domestic  waste  there  are  two  types  modeled, 
trickling  filter  and  activated  sludge.  ITYPE  supplies  the  program  with  this 
i nformation 

ITYPE  = 0 trickling  filter 

ITYPE  = 1 activated  sludge 

ITYPE  is  read  in  with  the  following  statements: 

READ  25,  ITYPE 
25  FORMAT,  ( 1 1 ) 

The  card  layout  for  ITYPE  =1  is 

1 3 3 V 5 

'mil// 

The  domestic  treatment  data  is  read  in  on  two  cards.  Appropriate 
numbers  are  entered  for  the  particular  type  of  treatment  used.  The  appro- 
priate data  for  each  type  is  as  follows: 

TRICKLING  FILTER 

AREP  primary  clarifier  area  in  acres 
TVOL  trickling  filter  volume  in  acre-feet 
TEMP  wastewater  temperature  in  °C 

R2  recycle  from  filter  effluent  (MGD) 

R3  recycle  from  secondary  clarifier  (MGD) 

Rl)  recycle  to  filter  (MGD) 

TFAREA  trickling  filter  area  in  acres 
K20  “ 0.23  for  1"  rack  media 


= 0.13  for  2V1  rock  media 
AREAS  secondary  clarifier  area  in  acres 
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ACTIVATED  SLUDGE 


ARE  P 

as  above 

AREA 

as  above 

ABVOL 

aeration 

basin  volume  (MG) 

MLSS 

mixed  liquor  suspended  solids  (mg/1) 

YN 

net  yield 

from  plant  log  data 

Rl 

act  i vated 

sludge  recycle  (MGD) 

The  data  are  read  in  with  the  following  statements: 

READ  26,  TEMP,  AREAP,  AREAS,  ABVOL , TFAREA,  TFVOL 

26  FORMAT  (F9.3,  IX,  F9.3,  IX,  F9-3,  IX,  F9 . 3 , IX,  F9.3.  IX,  F9.3) 

READ  27,  MLSS , YN , R1 , R2 , R3 , Rb 

27  FORMAT  (F9.3,  IX,  F9.3,  1 X , F9.3,  IX,  F9.3,  IX,  F9.3,  IX,  F9.3) 

The  card  layout  illustrating  the  first  two  of  these  values  TEMP  = 

20°C  and  AREAP  = 5 acres  is  as  follows: 

/ -a  3 V r <m  7 s f /£  ti  t% _____ 

I I I \l\o\-\o\o\o\  I II  1 1sM4|o|o| 

Next  the  flagging  data  discussed  above  is  read  in  by  the  FLAG(P) 
subroutine.  First  the  X flagging  level  is  read  in  for  each  parameter  of 
the  survey  in  the  same  parameter  order  used  in  previous  data.  All  X's  are 
read  in  on  the  same  card  with  the  statements: 

READ  (1,900)  (FLGPT  (I) , I - 1 , (NP) 

900  FORMAT  (25F3.2) 
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and  for  flag  levels  of  .10,  .10,  .10,  .99,  .05  etc.,  the  card  layout  is: 


The  manual  flagging  is  done  in  two  parts,  the  flags  for  the  sources 
and  the  flags  for  the  branches.  First  a card  reads  in  NF  the  number  of 
flags  for  parameters  to  be  measured  at  sources,  then  the  cards  which 
designate  the  source  number  N X ( I ) and  the  parameter  number  of  that 
source  NY(I)  as  determined  by  the  consecutive  ordering,  not  the  master 
list  number.  One  card  is  used  for  each  pair  to  be  designated.  The 
NF  is  the  number  of  pairs  to  be  read  in.  The  statement  and  the  card 
1 ayout  fo 1 low: 

READ  (1,  903)  NF 
903  FORMAT  (12) 

For  9 source  flags  the  card  layout  is: 

1 2 3 A 5 6 7 

DO  60  I = 1 , NF 
60  READ  (1 , 904)  NX ( I ) , NY  ( I ) 

90A  FORMAT  ( 12 , IX,  12) 

The  card  flagging  source  3 parameter  11  follows: 

12  3^567 


I 131  Mill  I I 


1 171  1 1 1 


D I 81  1=1,  NFB 

81  READ  904,  NX  (I)  , NY  ( I ) 


904  FORMAT  ( 12,  IX,  12) 

The  card  flagging  branch  5 parameter  12  follows 


1 2 3 4 5 6 7 
/ 1 151  11121  1 


This  completes  the  data  called  by  the  START  and  TOP  and 


the i r subrout i nes . 


61 


B 

RESOURCE  MATCHING  MODEL  USER'S  GUIDE 

1 . Data  Requ i red  : 

Data  required  by  the  resource  matching  model  can  be  divided 
into  two  classifications.  First,  there  is  data  generated  in  the 
process  and  topology  models;  and  secondly,  there  is  data  that 
must  be  supplied  by  the  user. 

The  data  generated  in  the  process  and  topology  model  are 
the  parameter  concentrations  or  other  measure  number  for  each 
parameter  for  each  branch  of  the  system.  These  values  are 
computed  for  a given  set  of  process  states  and  are  available 
to  the  resource  matching  program  and  its  subroutines  through  the 
common  dimension  statements. 

The  data  to  be  supplied  by  the  user  consists  of  the 
f o 1 lowing: 

*1.  The  number  of  measure  methods  available  for  each  parameter 
(NMA(N) ) 

*2.  The  sampling  frequency  for  each  parameter  (SAMFRE(N)) 

3.  The  number  of  pieces  of  equipment  in  all  resource  data 
(NPOFEQ) 

A.  The  names  of  the  pieces  of  equipment  and  time  constraints 
on  each  piece  (EQNAME) , (EQTIME) 

5.  The  data  defining  the  measurement  methods  for  each  para- 
meter (PMDATA) 

6.  Constraints  other  than  on  pieces  of  equipment 

7.  Sampler  assignments 

8.  Flow  measurement  data 


This  information  was  supplied  to  a subroutine  of  the  TOP 
program,  CHK1 , but  it  is  used  here. 


1 . EXPD  SCALE  PH  MET 


k. 


2.  VIS  SPECTROPHOT 

3.  PH  PREPARATION 

4.  PH  MEASUREMENT 

5.  AUTOMATIC  BURET 

6.  HEATER,  EXTRACT  RACK 

7.  CONDUCTIVITY  METER 

8.  ANAL  BALANCE 

9.  OVEN  (I05C) 

10.  Oven  ( I 80C) 

11.  DESSICATOR 

12.  TECHNICON 

13.  CARBON  ANALYZER 

14.  BLENDER 

15.  LAB  HOOD 

16.  KJDL  DIGEST  RACK 

17.  VACUUM  SOURCE 

18.  MAGNETIC  STIRRER 

19.  NEPHELOMFTER  (HACH) 

20.  SEPERATORY  FUNNELS 

21.  FUNNEL  RACK 

22.  MUFFLE  FURNACE 

23.  AA  SPECTROPHOT 

24.  TECH  AUTOANALYZER 

25.  GAS  CHROM  Wl  FID 

26.  CENTRIFUGE 

27.  HACH  TURBI DMETER 

28.  HOTPLATE 

EQUIPMENT  INCLUDED  IN  THE  RESOURCE  MODELS 
TABLE  III 


k. 
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For  each  parameter  to  be  measured  there  are  1 to  3 methods 
defined  for  making  the  measurement.  The  different  methods  may 
involve  different  techniques,  different  equipment  or  different 
personnel.  The  number  of  methods  available  for  each  parameter  is 
needed  and  is  read  in  as  NMA(N). 

A sampling  frequency,  SAMFRE(N),  is  required  for  each  parameter. 

This  is  the  number  of  samples  measured  per  dav  per  measure  point  for 
the  given  parameter.  When  the  commitment  of  equipment  and  personnel 
per  sample  is  combined  with  the  total  number  of  sample  points  and  the 
samples  per  day  for  the  length  of  the  survey  the  total  commitment  for 
that  parameter  is  determined  for  the  survey. 

The  resources  available  for  use  in  planning  a survey  are  defined 
in  a resource  model.  The  resource  model  defines  one  to  three  measure- 
ment methods  for  each  parameter.  Each  measure  method  model  lists 
the  following  information  about  that  method  in  matrix  form:  (Figure  IX) 

1.  Parameter  name 

2.  Method  name  and  number 

3.  Minimum  acceptable  concentration 

4.  Equipment  items  (name  and  number) 

5.  Time  commitment  of  personnel  and  equipment  for  sample 

6.  Space  commitment 

7.  Miscellaneous,  including  cost,  set-up  time,  location  of  tests, 
and  time  period  in  which  tests  must  be  run. 

This  data  is  supplied  to  the  program  in  the  Resource  Matching  subroutine. 
The  cards  defining  measure  methods  for  each  parameter  should  be 
considered  a library  of  methods  and  the  library  can  be  updated  by  including 
the  latest  or  best  values  on  the  method  cards  used. 


PARAMETER  MEASURE  METHOD  DATA 
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The  program  reads  resource  information  in  as  follows: 

1.  The  number  of  pieces  of  equipment  included  in  all  measure 
methods  is  read  in.  (NPOFEQ) 

2.  The  names  of  the  pieces  of  equipment  are  read  in  along  with 
the  total  amount  of  time  available  for  the  survey  on  that 
piece  of  equipment  (the  order  in  which  pieces  of  equipment 
are  read  in  assigns  them  a number).  EQNAME(I),  EQT I ME ( I ) . 

3.  The  parameter  method  data  is  now  read  in  in  matrix  form. 

One  set  of  five  cards  for  each  method  for  each  parameter. 

The  first  card  for  each  parameter  gives  the  method  name  and 
the  minimum  concentration.  (MENAME)  (PCRM) . The  next 
four  cards  read  in  the  five  numbers  on  each  line  of  the 
PMDATA  model  illustrated  in  Figure  IX.  One  line  is  read 

in  per  card.  (PMDATA) 

The  next  resource  matching  data  required  is  a listing  of  con- 
straints other  than  those  already  listed  for  each  piece  of 
equipment.  These  constraints  are  given  on  a total  survey  basis 
and  are:  Total  Van  Space,  analysts  (1)  time,  analyst  (2)  time, 

analyst  (3)  time,  analyst  (4)  time  and  cost.  The  different 
classes  of  analysts  have  been  defined  by  AEHA  Chemists  in  terms 
of  abilities  and  preferences  for  performing  certain  functions. 
Analyst  (Class  I),  for  example,  would  be  used  to  collect  samples 
and  subdivide  them  into  test  size  samples.  Analyst  (Class  2) 
could  perform  most  routine  analyses.  Analyst  (Class  3)  could 
perform  the  more  complicated  analyses  and  would  be  used  on 
routine  analyses  only  as  second  choice.  Analyst  (Class  4)  mostly 
supervisory  but  could  fill  in  for  Class  2 and  3 if  needed. 

If  the  data  for  the  van  space  is  properly  defined,  this 
information  can  be  used  to  determine  if  additional  facilities 


^4 


will  be  needed  for  the  survey,  (e.g.  , additional  wet  chemistry 
laboratory  facilities). 
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Other  resources  in  the  survey  separate  from  the  analytical 
chemistry  are  the  resources  required  for  collecting  the  samples. 
These  resources  include  sampler  equipment,  equipment  time,  and 
personnel  time  commitment  to  collect  the  samples.  The  information 


requ i red 

by  the 

SAMPLE  subroutine 

is  shown  in  Table  IV. 

The 

personne 1 

time 

requires  a commitment  of  analyst  (Class 

1 ) time. 

tt 

Sampler  Name 

Time/Sample 

Set-Up  Time 

1 

EASY  GRAB 

5.0 

5.0 

2 

HARD  GRAB 

15.0 

5-0 

3 

ISCO 

20.0 

5-0 

A 

PROTECH 

5.0 

5-0 

SAMPLER  DATA 
TABLE  IV 

In  addition  to  the  sampler  data  of  Table  IV,  the  survey 
officer  must  decide  during  the  preliminary  survey  which  samplers 
will  be  used  at  each  branch  if  the  survey  plan  calls  for  a 
measurement  at  that  branch.  This  data  is  supplied  to  the  SAMPLE 
subroutine  by  reading  in  the  appropriate  sampler  number  from  the 
left-hand  column  of  Table  IV  for  each  branch.  NROUT(l). 

The  last  items  of  data  required  are  for  the  flow  measurement. 
The  flow  measure  method,  time  per  measurement  and  set-up  time  are 
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# 

1 

2 


Flow  Meter  Time/Meas. 

Stevens  Recorder 
Gurley  Meter 


Set-up  Time 
20.0 
30.0 


FLOW  METER  DATA 
Table  V 

This  information  is  supplied  to  the  SAMPLE  subroutine.  In 
addition,  the  survey  officer  must  decide  for  each  branch  if  a 
flow  measurement  is  to  be  made  and  if  the  measurement  is  to 
be  made  will  meter  type  #1  or  #2  be  used.  This  information  is 
supplied  to  the  program  by  reading  in  for  each  branch  a zero 
for  no  measurement,  a one  for  a measurement  using  a Stevens 
recorder,  or  a two  for  a measurement  with  a Gurley  meter. 

I I . Data  Format 

The  formats  and  card  layout  for  reading  in  the  resource 
matching  data  is  given  next. 

The  first  card  read  in  by  the  resource  matching  model  is 
the  one  giving  the  number  of  pieces  of  equipment  in  PMDATA. 
This  is  NPOFEQ  and  is  read  in  by  statements 
READ  135,  NPOFEQ 
135  FORMAT  (13) 

the  card  layout  for  32  pieces  is 
12  3^5 

a 13  i2i  i i r 


The  next  group  of  cards  list  equipment  names  and  the  time 


* 
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available  on  that  item  of  equipment.  Two  items  of  equipment 
per  card  (one  item  shown).  See  Figure  XIX. 

READ  145.  ( (EQNAME  (I . J)  J - 1,5),  EQTIME  (I),  I = I,  NPOFEQ) 

145  FORMAT  (2(5A4,  F10. 1 , 1 OX) ) 

1234567  20  30  40 


| E | X | P | D | | S | C | A | L | E | | P | H | |M|E|T| 


5 | 8 | 8 | 0 | . J 0 


The  next  set  of  data  cards  contain  the  parameter  measure- 
ment method  data  (PMDATA) . 

Five  cards  are  used  for  each  parameter  for  each  method. 

The  first  in  each  group  of  five  gives  the  method  name  and  the 
minimum  concentrations  required  by  the  method.  The  read  and 
format  statement  are: 

READ  410,  (MENAME  (N , J,  II),  II  = 1,5),  PCRM  (N,j) 

410  FORMAT  (5X,  5A4,  F 1 0 . 5) 

Example  Card 

5 25 

/ 

| | | | |E|L|E|C|T|R|0|M|E|T|R|l|C|  ||  | | | | | | | | | . | 0 ; 1 0 0 


MENAME  ( N , J ) 

Nth  parameter  Jth  method 


PCRM  ( X , J ) 


The  next  four  cards  supply  the  four  lines  of  data  of  the 


parameter  measure  method  matrix;  One  line/card  in  proper 
order.  The  read-in  and  format  statement  are: 


0030  L = I ,4 

30  READ  430,  (PMDATA  (N,  J,  L,  LI),  LI  - 1,5) 
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430  FORMAT  (F5.0,  F10.2,  Fl 0. 3 , F10.2,  F 10. 3) 

Example  card  for  first  line  of  data 

15  15  25  35  45 

A I I 4 | • | 0 | | 1 | . | 0 | 0 | +-n-p  j 3 | . | 0 | 0 | 0 | | 0 J . | 0 | | 0 | . | 1 | 0 | 

EQ  * T/SAMP  VAN/SP  ANAL/TIME  COST 

The  other  three  lines  are  read  in  on  the  next  three  cards 
as  L iterates  from  1 to  4. 

If  there  is  a second  method  (J)  for  the  same  parameter  (N), 

PMDATA  is  read  in  for  the  next  method  (five  cards).  The  next 
parameter  is  then  considered  as  N goes  to  N + 1. 

After  all  of  the  PMDATA  cards,  the  constraint  card  follows. 

This  lists  any  constraints  on  Van  space,  analyst's  time,  costs  or 
other.  This  program  is  set  up  to  crnsider  six  constraints  in 
addition  to  the  constraints  on  equipment  time. 

The  read  and  format  statement  are: 

READ  650,  (CNSTAR  (M) , M » 1 ,6) 

650  FORMAT  (4F  20.3) 

Example  card 

20  40  60  80 

/ 1 1 | 0 I 0 I 0 I . I 0 I 0 I 0 |*2,|6|3|  0 I 0 I . |0|0|0|***|6|  3|0|0|  . |0|0|0|^  1 1 2 | 6 [ 0 . 

VAN  SPACE  ANAL  1 ANAL  2 ANAL  3 


A second  card  using  the  same  format  is  required  to  list  ANAL  A 
and  cost  constraints. 


The  format  for  reading  in  the  sampler  name,  sample  time 
and  set-up  time  is  illustrated  by  the  statements 
D030  1-1,4 

30  READ  903,  (NSAMPL  (I,  K) , K - 1,5)  SAMPTI  (I),  SET  UP  (I) 
903  FORMAT  (5AA,  2F6.  1) 


An  example  card  for  reading  in  EASY  GRAB,  5.0  minutes  per 
sample,  and  5.0  minutes  set-up  time  is  shown 


1234S6789 


26 


32 


/ 

| E | A | S | Y | | G | R | A | B | 


5 ] • | 0 | 


51-10 


The  listing  of  the  sampler  type  to  be  used  for  each  branch 
If  a sample  is  to  be  taken  at  that  branch  is  read  in  by  the 
following  statements 

READ  906,  (NROUT  (I) , I - 1 ,NB) 

906  FORMAT  (4012) 

The  card  layout  for  reading  in  sampler  1 for  branch  1, 
sampler  3 for  branch  2,  etc. 

1 2 3 4 5 6 7 8 

I 111  1 3 1 |1|  | 2 | 


The  format  for  reading  in  the  flow  meter  name,  sample  time, 
and  set-up  time  is  given  by  the  statements  (this  is  an  extension 
of  the  sampler  read  in). 
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time  and  20.0  minutes  set-up  time  is  shown 

123456  26  32 


The  listing  of  the  flow  measure  method  to  be  used  for  each 
branch 

READ  906,  (NFLOW(I) , I = 1 , NB) 

906  FORMAT  (4012) 

Example  card  with  no  flow  reading  on  branches  I,  2,  A,  5,  and 
a Stevens  recorder  on  3 and  a Gurley  meter  on  6. 
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MODEL  VERIFICATION  AND/OR  UPDATING 


I . INTRODUCTION 

The  model  verification  and  or  updating  is  done  using  the  CORRCT  subroutine. 
The  CORRCT  subroutine  accepts  site  measurement  data,  compares  it  with  model 
predicted  measurement  data  and  computes  and  prints  out  suggested  new  values 
for  the  parameters  of  the  process  models. 

The  subroutine  CORRCT  is  used  in  conjunction  with  the  program  sub- 
routines which  make  up  the  site  model.  They  are  MAIN,  PROCES,  START,  TOP, 
NEUTRA,  SETTLE,  DOME,  CHKl , and  CHK  2. 

The  site  model  package  with  CORRCT  has  two  objectives.  First,  the 
process  models  can  be  updated  or  improved  using  preliminary  survey  data. 

While  the  best  data  available  were  used  in  constructing  the  process  models, 
there  was  considerable  variation  in  data  for  a given  process  from  one  plant 
to  another  or  from  one  season  to  another.  Better  process  model  values  may  be 
obtained  with  CORRCT  using  measurements  from  judiciously  selected  grab 
samples  taken  on  a preliminary  survey.  The  output  of  CORRCT  lists  new  values 
for  certain  process  parameters.  If  the  survey  officer  wishes,  he  may  enter 
these  new  values  into  the  process  models  using  procedures  outlined  in  the  dis- 
cussion of  the  PROCES  subroutine. 

A second  use  of  the  site  model  package  with  CORRCT  would  be  verification 
of  the  model.  After  the  survey,  measurements  are  available  for  obtaining 
max,  min,  average  or  other  expected  values  for  each  parameter  at  each  branch 
for  known  levels  of  production.  These  values  can  be  entered  into  the  site 
model  package  with  CORRCT  to  determine  values  for  process  model  parameters 
which  can  be  used  to  verify  the  model  or  refine  the  model  as  appropriate. 
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With  confidence  in  the  model,  one  can  now  use  it  to  predict  levels  of  potential 
pollutants  for  production  levels  other  than  those  which  occurred  during  the 
survey.  This  should  give  added  significance  to  the  survey  report. 

Ideally,  if  measurements  are  to  be  made  to  verify  the  process  models, 
the  measurepo i nts  should  be  chosen  such  that  each  process  will  be  determined 
uniquely  by  the  measurements.  However,  if  the  survey  officer  who  determines 
the  measure  points  does  not  visualize  the  system  in  terms  of  the  same  processes 
as  are  used  in  the  process  model,  measurements  may  not  uniquely  determine 
process  parameters.  The  program  is  designed  to  accept  such  non-unique  mea- 
surement data  and  uniquely  determine  processes,  where  possible,  and  make 
approximations  on  other  processes  informing  the  user  of  the  determination  is 
non-un i que . 

The  measure  points  for  uniquely  determining  all  processes  can  be  identi- 
fied by  referring  to  the  linear  graph  of  the  system  (Figure  III  for  example) 
which  display  all  processes  and  all  possible  measure  points.  On  the  linear 
graph,  identify  all  of  the  vertices  (junctions)  nearest  to  the  processes. 

If  at  a vertex  there  are  N branches  incident  to  a vertex,  then  measurements 
must  be  taken  or  implied  at  N- 1 of  these  branches.  By  taking  N-l  measurements, 
the  N*"  measurement  is  implied.  Note  that  one  of  the  N-l  measurements  at  a 
vertex  A may  be  implied  by  measurements  at  a vertex  B or  other  vertices. 

I I . INPUT  DATA  AND  OUTPUT  DATA 

After  the  decisions  are  made  which  determine  where  measurements  are  to 
be  made,  this  information  is  supplied  to  the  CORRCT  subroutine  in  the  form 
of  the  MESUR(I)  matrix.  The  MESUR(I)  matrix  contains  a "0"  for  each  branch 
where  measurements  are  not  taken  and  a "1"  for  each  branch  where  measurements 


are  taken. 


It  may  not  be  appropriate  or  even  possible  to  check  every  parameter  in 
a process  model.  The  subroutine  CORRCT  only  applies  to  conservative  parameters 
[i.e.  only  those  measured  in  mass  terms  (mg/1)  or  flow].  The  program  allows 
the  user  to  select  the  one,  two,  or  several  (up  to  20)  parameters  to  be 
included  in  each  model  verification.  The  parameters  to  be  considered  in  a 
particular  study  are  identified  by  their  numbers  in  the  master  list  of 
Table  II.  The  Flows  must  be  included  in  any  branch  where  other  parameters  are 
cons i dered . 

The  number  of  parameters  and  the  parameter  numbers  are  supplied  to  CORRCT 
for  each  branch  where  measurements  are  to  be  taken  using  BPLIST  ( I B , N ) , where 
IB  is  the  branch  number  and  N is  the  parameter  number.  The  first  number  N 
for  each  branch  is  the  number  of  parameters  considered  in  that  branch.  The 
subsequent  numbers  N are  the  parameter  numbers  for  that  branch. 

The  next  inputs  to  CORRCT  are  the  measured  parameter  values  YM(|B,N). 

For  each  branch  IB  the  first  value,  for  N=l , must  be  the  flow  in  MGD.  The 
remaining  YM(|B,N)  input  for  a branch  are  the  measure  values  of  the  parameters 
listed  in  BPLIST(IB.N). 

The  subroutine  CORRCT  output  lists  each  source  (process)  by  number, 
notes  if  the  determination  is  unique  or  estimated,  and  gives  the  old  parareter 
and  the  suggested  new  parameter  for  each  parameter  verified.  If  the  deter- 
mination is  unique,  the  branch  where  the  defining  measurement  was  taken  is  listed. 

III.  DATA  FORMAT 

The  data  required  by  CORRCT  is: 

1.  MESUR(I)  , indicating  if  a measurement  is  made  in  branch  I 

2.  BPLIST(lB.N)  listing  for  N=l,  the  number  of  parameters  to  be 


verified  in  this  branch  and  for  other  N is  the  parameter  numbers 
for  parameters  to  be  verified  for  branch  IB. 


75 


3.  YM(IB,N),  listing  for  N=1 , the  flow  in  branch  IB  and  for  other  N's 
the  measurement  values  for  each  parameter  listed  in  BPLIST  ( I B,N) . 
MESUR(I)  is  read  in  by: 

READ (1,1)  (MESUR(I),  1=1,  NB) 

1 FORMAT  (80  I 1) 

Card  layout  for  measurements  in  branches. 

1,  2,  3,  5,  7,  9. 

/ J 3 q JT  6 7 / 0 

/ I'lUUI  111  U\  111  I 

BPLIST(IB,N)  and  YM(lB,N):  are  read  in  one  branch  at  a time  by  the  state- 
ments. 

DO  2 IB  = 1 , NB 

IF  (MESUR(IB) . EQ.  0)  GO  TO  2 

READ  (1,3)  (BPLIST  (IB,N),  N = 1,20 

3 FORMAT  (2012) 

I BPLIST  = BPLIST  ( I B , 1 ) 

READ  (1,4)  YM  (I B , N ) , N = 1,  BPLIST 

4 FORMAT  (F  20.6) 

Card  layout  for  two  measurements  on  a branch  IB,  one  flow  and  the  other 
acidity.  Note  the  parameter  number  for  flow  does  not  exist  and  is  not  given 
in  BPLIST.  (Acidity  number  is  2k) 

BPLIST  Card 

' 1 3 v T 6 

[MIlhlYl  I 1 


YM  card  for  flow  = 0.78  and  acidity 
/ -2  3 V 5 '6  7 g 1 /O 

(i  M'l’UI  I II  I l 


12  3 V s'  b 1 5 •/  /C 


I I Uls'l'M 
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Review  of  Input  Data 

The  Diagrams  on  the  next  three  pages  serve  as  a review  of  the 
data  required  by  the  Survey  Planning  Program  and  also  it  serves  to 
emphasize  that  the  order  in  which  the  data  is  read  in  must  be  exactly 
as  shown  in  these  diagrams. 
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ORDERING  OF  INPUT  DATA 


MASTER  = largest  parameter  number 


NFLG  Flag  to  bypass  RM 
EPSLON  Flag  to  call  CORRCT 


N the  number  of  procener  in  the  plast  model 


XNAME  (JM) , CAP  (JM)  process  model  name  and  capacity  and 
card  per  process  in  topological 


NOD  Flag  to  indicate  modifications 


NMOD  The  number  of  modifications 


YNAME  (NM) , EI.E(MN),  VAL  (NM)  Process  name,  parameter 

member,  parameter  value.  One  card  per  modification 


NPRO  Flag  to  indicate  complete  processes  to  be  changes 


NUMP  Number  of  parameter  in  the  master  list 


ZNAN1E  (JN)  , FEW  (JN)  Name  and  Flo\.  of  processes  to  be  changed 


SPLIST  (JF,,  JK)  Parameter  value,  one  card  for  each  parameter 
in  the  master  list  for  each  changed  process. 


NS,  NR,  NP,  NTOP  The  number  of  sources,  branches,  parameter-, 
and  outfalls 


NFLIST  (I)  List  of  the  pol Intent  number  used 


A (J,I)  The  topology  matrix,  one  card  per  rov 


NBRANCH  branch  names 


LENGTH  Length  of  survey  is 


ISTOP  The  nubmer  of  nor -compe ting  parameters 


POLN,  NMA,  SAMFRE  The  names  of  all  parameters  in  the  master 
up  to  MASTER,  the  number  of  measure  methods  avoilable 
and  the  sampling  frequency 


Treatment  Facilities  Data 

*KEA,riNEUT  ,MEAE  ,UCMS  ,MACY  , I TYPE  ,TEnP  ,ni\E*P  , AREAS  ,«oV0L  ,TFAKEA  ,TFVOL  , 
mlss!ri.r2.’r3.r4. ! 

FLGPT  The  flagging  level  is  read  in  for  each  parameter 


The  number  of  manual  source  flags 


NX,  NY 


The  source  and  parameter  flagged  one  card  for  each  \'F 


The  number  of  manual  branch  flags 


NX,  NY  The  branch  and  parameter  flagged  one  card  for  each 
NFB 


NPOFEQ 


The  number  of  pieces  of  equipment  in  the  resource  data 


EQNAME , EQTIME 
survey 


Equipment  name  and  time  avoilable  during 


MENAME , PCRM  Method  name  and  minimum  concentration  for  each, 
parameter  for  each  method 


PMDATA  Method  resource  data 
for  each  method 


- four  cards  for  each  parameter 


CNSTAR  Lists  constraints  on  van  space,  analysts  time, 
two  cards 


etc  , 


NSAMPLE,  SAMPTI , SETUP  Sampler  name,  time  per  sample, 
up  time  for  each  type  of  sampler 


set 


NROUT  List  the  sampler  to  be  used  for  each  branch  if  a 
sample  may  be  taken  there 
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X SAM PL , SAMPTI , SETUP  Flow  meter  name,  sample  time  and 
set  up  time  for  each  type  of  meter 


X F L Oh  Lists  the  flow  meter  to  be  used  for  each  branch 
where  a reading  may  be  taken 


MEASUR  Lists  branches  where  measurement  data  is  taken 


BPLIST  For  each  branch  measured,  give  the  number  of  parameters  to  be 
verified  and  lists  them  by  numbers. 


YM  Gives  the  appropriate  measured  values. 


Output  Data 


The  output  data  falls  into  two  general  c 1 ass i f i cat  ions  , 
output  which  affirms  the  input  data  and  output  which  is  a 
result  of  analysis.  In  this  section  the  output  data  will  be 
reviewed.  If  the  data  is  a reiteration  of  input,  it  will  be 
identified  and  no  further  discussion  will  be  given.  Other 
output  will  be  considered  in  more  detail. 

All  of  the  output  data  shown  on  Figure  X are  reiterations 
of  input  data.  The  first  10  lines  give  the  mnemonics  of  the 
processes  in  the  survey  and  their  operating  levels.  The  "1" 
in  line  eleven  indicates  that  a process  model  is  to  be  altered. 
The  "1"  in  line  12  indicates  that  there  is  to  be  one  altera- 
tion. Line  13  defines  that  change.  The  "1"  in  line  1A  indi- 
cates that  there  are  new  processes  to  be  read  in. 

The  output  data  shown  on  Figure  XI  are  mostly  reiterations 
of  input  data.  The  first  line  is  NP , the  number  of  parameters 
to  be  measured  in  the  survey.  The  second  line  indicates  that 
the  survey  planning  period  is  one  day.  The  remaining  17  lines 
list  the  parameters  to  be  measured  on  this  survey,  their 
assigned  number,  (left  hand  column),  NMA  and  SAMFRE.  The 
last  column  on  the  right  gives  the  total  number  of  samples  per 
point.  This  is  the  product  of  sample  frequency  and  the  length 
of  the  survey  in  days. 

Figure  XII  shows  the  next  groups  of  output  data.  All  of 
these  are  reiterations  of  input  with  an  ordering  assignment. 
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for  the  processes.  At  the  top  of  the  page  three  numbers,  NS, 
NB,  and  NTOP  are  given.  The  next  block  of  output  is  a re- 
listing of  the  topology  matrix,  A(j,l).  The  last  group  in 
this  figure  is  headed  pollutant  matrix  data  which  lists  the 
process  names,  their  assigned  number  and  their  flows  in  MGD. 

The  next  block  of  data  is  illustrated  in  Figure  X 1 1 1 . 

It  is  a listing  of  each  source  mnemonic,  its  assigned  number 
and  the  parameter  value  for  each  of  the  pollutants  for  each 
process.  This  is  essentially  a reiteration  of  the  input  pro- 
cess data. 

The  output  data  in  Figure  XIV  illustrates  the  flagging 
information.  The  first  block  of  data  gives  the  flagging 
levels  for  each  parameter  for  automatic  flagging.  The  second 
block  of  data  is  a list  of  the  manual  flags  at  sources.  Not 
shown  but  also  part  of  the  output  is  a listing  of  the  manual 
flags  at  the  branches. 

Figure  XV  lists  the  assigned  branch  numbers  and  names 
and  also  includes  the  calculated  flow  for  each  branch. 

Figure  XVI  lists,  at  the  top  oc  the  page,  the  number  of 
branches,  NB,  and  the  number  of  parameters,  NP,  in  the  study. 
The  remaining  data  on  the  Figure  XVI  illustrates  the  results 
of  the  TOP  calculations.  For  each  branch  the  name  and  number 
of  the  branch  is  listed  and  along  with  this  the  parameter  value 
is  listed  for  each  pollutant  in  the  study.  The  negative  signs 
on  certain  of  the  numbers  is  not  to  be  interpreted  as  indi- 
cating negative  values.  Rather,  it  is  an  indication  that 
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this  parameter  is  flagged  for  that  particular  branch.  Figure 
XVI  gives  only  an  illustrative  sample  of  this  output  data. 

Figure  XVII  lists  the  mass  in  pounds  per  day  that  passes 
through  each  branch  of  the  system.  This  mass  number  is 
meaningful  for  such  parameters  as  TOC  and  NITRITE/NITRATE  but 
has  no  meaning  in  the  case  of  pH  and  Conductivity.  Again  the 
negative  numbers  are  flagged  points  and  not  negative  quantities. 

Figure  XV III  shows  the  ordering  of  the  measure  points  as 
done  by  the  level  subroutine.  This  data  tells  the  number  of 
measure  points  required  to  detect  all  of  a parameter  if  meas- 
urements are  made  at  levels  different  from  the  outfall.  For 
example,  if  the  measurements  are  to  be  made  at  level  2 then 
there  are  three  measure  points  required. 

Figure  XIX  is  a reiteration  of  input  data  including 
NPOFEQ,  equipment  name  and  number,  and  time  available  on  each 
piece  of  equipment. 

Figure  XX  illustrates  the  next  groups  of  data  read  out. 

This  is  the  PMDATA,  the  measurement  method  data  for  each 
measure  method.  This  is  read  out  for  each  parameter  for  each 
method . 

Figure  XXI  is  a reiteration  of  the  constraint  data  input. 

The  next  output  data  is  a listing  of  each  branch  of  the 
system  and  the  parameters  that  are  to  be  measured  at  each 
branch.  Figure  XXII  illustrates  this  output  for  several 
branches . 

Figure  XXII I lists  the  sampler  data  at  the  top.  The 
figure  also  lists  computed  sampler  commitment  for  each  branch 
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as  a result  of  the  completed  survey  plan. 

Figure  XXIV  lists  the  flow  meter  data  at  the  top  of  the 
page.  The  figure  also  lists  computed  flow  meter  time  commit- 
ments as  a result  of  the  completed  survey  plan. 

The  next  group  of  data  is  illustrated  by  Figure  XXV. 

This  output  lists  the  following  data  for  each  parameter 
measured : 

(1)  The  parameter  name  at  the  upper  left  corner. 

(2)  At  the  upper  right  corner  the  output  lists  the 
branches  in  which  the  measurements  are  to  be 
made  in  the  survey,  the  measure  method  to  be 
used  at  each  branch,  and  the  expected  value 

of  the  parameter  at  that  branch.  The  asterisks 
indicate  flagged  measure  points. 

(3)  The  total  number  of  samples  of  that  parameter 
during  the  survey  planning  period. 

(4)  The  items  of  equipment  required  in  measuring 
this  parameter,  the  total  time  commitment  on 
each  for  these  measurements,  and  an  indication 

if  a time  constraint  was  violated  (for  the  survey, 
as  a whole, does  not  indicate  at  which  parameter 
the  violation  occurred). 

(5)  The  time  commitment  of  each  of  the  types  of  analyst 
in  making  these  measurements  and  an  indication  of 


time  constraints  violated. 


Figure  XXVI  shows  the  read  out  data  which  is  a summing  up 
of  total  commitments  of  all  the  equipment  and  analysts  for 
the  total  survey  planning  period.  Along  with  the  total  com- 
mitments, violations  are  indicated  and  the  amount  noted. 

Figure  XXVII  shows  the  output  data  of  the  CORRCT  subroutines. 
When  CORRCT  is  used  the  usual  site  model  output  from  TOP  is  read 
out  and  in  addition  the  data  in  the  form  of  Figure  XX V | | is  output 
to  give  the  suggested  new  values  for  process  parameters.  There 
are  two  types  of  outputs.  In  one  type  where  the  values  are 
determined  directly  from  measurements,  values  are  listed  as  unique 
determinations.  For  other  types,  where  the  values  cannot  be 
determined  directly  from  measurements,  values  are  listed  as 


estimated  values. 
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SOURCE  NUMBER  I IS  UNIQUELY  DETERMINED  ON  BRANCH  1 WHERE  3 PARAMETERS  WERE  MEASURED. 
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